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ABSTRACT 

We present the results of a deep X-ray survey conducted with XMM-Newton, centred on the 
UK ROS AT 1 3 H deep field area. This region covers 0.18 deg 2 and is the first of two areas cov- 
ered with XMM-Newton as part of an extensive multi-wavelength survey designed to study the 
nature and evolution of the faint X-ray source population. We have produced detailed Monte- 
Carlo simulations to obtain a quantitative characterisation of the source detection procedure 
and to assess the reliability of the resultant sourcelist. We use the simulations to establish a 
likelihood threshold above which we expect less than 7 (3%) of our sources to be spurious. 
We present the final catalogue of 225 sources. Within the central 9', 68 per cent of source po- 
sitions are accurate to 2" making optical follow-up relatively straightforward. We construct 
the N(> S) relation in four energy bands: 0.2-0.5 keV, 0.5-2 keV, 2-5 keV and 5-10 keV. In 
all but our highest energy band we find that the source counts can be represented by a double 
powerlaw with a bright end slope consistent with the Euclidean case and a break around 10~ 14 
erg cirT 2 s . Below this flux the counts exhibit a flattening. Our source counts reach densities 
of 700, 1300, 900 and 300 deg~ 2 at fluxes of 4.1 x 10~ 16 , 4.5 x 10~ 16 , 1.1 x 10~ 15 and 5.3 x 10~ 15 
erg cirT 2 s in the 0.2-0.5, 0.5-2, 2-5 and 5-10 keV energy bands respectively. We have com- 
pared our source counts with those in the two Chandra deep fields and Lockman hole and find 
our source counts to be amongst the highest of these fields in all energy bands. We resolve 
> 51% (> 50%) of the X-ray background emission in the 1-2 keV (2-5 keV) energy bands. 
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1 INTRODUCTION 

It is widely accepted that the majority of the Cosmic X-ray Back- 
ground (XRB) arises from the integrated emission of discrete 
extragalactic sources ISchwartz et allll976tlGiacconi & Zamoranil 
Il987t lMaccacaroetalJll99ll) . The energy density of the XRB 
peaks at ~30 keV, but the first imaging surveys were carried 
out at much lower energies: <3.5 keV with Einstein and <2 
keV with ROSAT . By the lat e 1990s ROSAT surveys had re- 
solve d 70-80% of the soft XRB. iShanks et alJl99lllHasinger etai] 
I1993L 1 1998k iMcHardv et all Il998l) . Subsequently, deep XMM- 
Newton and Chandra surv eys have essentially resolved the soft 
XRB into discrete sources iMushotzkv et alj200(i Easinger et alJ 
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| 200J iBrandt et alJ 120011: iTozzi et alJ EoOlt iRosati et alJ l2002t 

lAlexander et alj |20031. Optical follow up of these sources has re- 
vealed a population composed primarily of unobscured broad line 
active galactic nuclei (AGN) , with an increasing fraction of ab- 
sorbed AGN at fainter fluxes iMcHardv et al l 1998 : ISdmiidtetalJ 
1 1998 } : IZamorani et aljfl999l Lehmann et aljboOlf : bzokolv etall 
l20031:lBarger et all2003h 

In order to investigate further the nature of the ob- 
scured population one has to conduct surveys at harder 
energies (>2 keV), which are l ess sensitive to abso rption. 
Surveys carried out using ASCA ^GeorgatrtOTOulo^UtlJ ^97 ; 
Cagnoni. Delia Ceca & Maccacarol Il998t lUedaetalJ 199fc . 
200 it llshisaki et alJ [200 ll) and Bepp oSAX iFiore et alJ ll99S 1 
Giommi et alJ l200d : iFiore et alJ l200lh resolved 25-35% of the 
XRB above 2 keV. More recently, the XMM-iVwto« and Chan- 
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dra deep field surveys ha ve resolved 60-90% of the hard (>2 
keV) X-ray background ( Hasin eer et all 1200 it I Giacco ni et alJ 
| 200lt iTozzi et al] 1200 it ICowie et alJ 120021: iRosati et alJ I2OO2T 
lAlexande^tiril2()o"3TlM^ fluxes a factor 

lOOx fainter than the ASCA and BeppoSAX surveys. The wide 
range in resolved fraction arises not only due to the variation in 
source counts between the surveys, but also due to the uncertainty 
of the absolute normalisation of the hard XRB, for exa mple the 
BeppoSAX XRB normalisation f romlVecchi et akHl999h is -30% 
higher than the ASCA value from Gendreau et al. ( 1995). 

Optical follow up studies of the deepest surveys find a pre- 
dominanc e of z < 1 objects which do not show broad emis- 
sion lin es 1 Barg er et alj200lHTozzi et alj200lHBarger et alj2002b 
lLamer et al 120031 : iBarger et all2003l) . There is an increasing con- 
tribution from normal galaxies at the faintest fluxes, and it ap- 
pears likely that normal galaxies will outnumber AGN below 0.5- 
2 keV fluxes of 10~ 17 erg cnr 2 s" 1 <Bauer et alJl2004l) . A signif- 
icant fraction of the hard X-ray sources in these fields are opti- 
cally faint, with R > 24 and are therefore difficult to identify opti- 
cally jAjexande^yii]200ll |McH^^^]20o3;|Aiexmd^^] 
120031) . 

Despite the great advances made in detecting increasingly 
fainter X-ray sources, the physical nature and evolution of the faint 
hard X-ray population remains largely unknown. The redshift dis- 
tribution and column density distribution of the absorbed AGN are 
still poorly constrained. 

Another issue which needs addressing is the relationship be- 
tween gas and dust absorption in AGN. There have been sev- 
eral cases of a mismatch between optical and X-ray classifica- 
tions, in dicating a wide rane e in dust/gas ratios for obscured 
sources ( Akivama et al. 2000; P age et al. ]200ll Comas tri et al] 
2221; iMaiolino et al] 12001 JlLoaring et al] 



2004). 



2003; Carreraetal. 



In particular, high quality X-ray spectra are needed to deter- 
mine the dominant X-ray emission mechanisms and the amount of 
absorption. We have therefore used XMM-Newton to carry out deep 
surveys of two widely separated fields to probe the X-ray popula- 
tion down to fluxes ~ 10~ 15 erg cnr 2 s in the 0.5-2 and 2-5 keV 
energy bands. The source count s in these e nergy bands exhibit a 
break at ~ 10 erg cm 2 s- 1 jRosati et alj|2002l) around which 
the maximum contribution to the XRB per logarithmic flux interval 
occurs. 

This paper presents the X-ray catalogue derived from the 
first of these two XMM-Newton surveys, carried out in the UK 
ROSAT deep field area (hereafter the 13 H deep field). The field 
has also been observed with a mosaic of Chandra pointings which 
cover the whole XMM-Newton field of view and provide accurate 
source positions. It is complemented with multiwavelength follow 
up in the UV, optical, near-IR, mid-IR and radio ISevmour et alJ 
2004). The Chandra c atalogue has already been presented else- 
where (McHardv et al. 2003); here we present the XMM-Newton 
catalogue and observed source counts. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 XMM observations 

The 13 H deep field is centred on the sky co-ordinates RA 13h 34m 
37.1s, Dec +37° 53' 02.2" (J2000). The XMM-Newton observa- 
tions were carried out in three separate revolutions during June 
2001 for a total exposure time of 200 ks. The European Photon 



Imaging Cameras (EPICs lTurner et all200ll) were operated in stan- 
dard full-frame mode. The thin filter was used for the PN camera, 
while the thin and medium filters were alternated for the MOS 1 and 
MOS2 cameras. TableQ gives a summary of the observations. 

The data were processed using the XMM-Newton Standard 
Analysis System (SAS) version 6.0. Approximately 40% of the 
total observation time was affected by high particle background 
flares, arising from soft protons hitting the detector. The data were 
therefore temporally filtered to remove these high background pe- 
riods. In practice, times where the 5-10 keV count rates exceeded 2 
s -1 in the MOS cameras and 4 s _I in the PN camera were excluded. 
Filtering reduced the total useful exposure time from ~200 ks to 
~120 ks. The net live times for the individual detectors after the 
periods of high background were excised are listed in TableQ 

A significant component of the EPIC background comes from 
instrumental emission lines, in particular the Cu Ka line at 8.1 
keV in the PN and the Al Ka line at 1.5 keV in both detectors 
jLumb et all2002l) . Events with energies close to those of the emis- 
sion lines were filtered out to minimise the instrumental contribu- 
tion to the background. Events in bad columns, bad pixels and ad- 
jacent to chip edges were also excluded. 

Images and exposure maps were then constructed from each 
observation for each detector in four energy bands: 0.2-0.5 keV, 
0.5-2 keV, 2-5 keV and 5-10 keV. Single-pixel events were used 
to construct the PN 0.2-0.5 keV image. Single, double and triple 
events were used to construct the higher energy PN images. For 
MOS, all valid event patterns were used to construct the images re- 
gardless of energy band. In each energy band, the exposure maps 
were scaled to the PN thin filter response. The images and ex- 
posure maps from the different detectors and observations were 
then summed to produce one image and one exposure map per en- 
ergy band. The response-weighted summation over each observa- 
tion and telescope gives total on-axis PN-equivalent live exposure 
times of 152 ks, 161 ks, 179 ks, and 160 ks, in the 0.2-0.5 keV, 
0.5-2 keV, 2-5 keV, and 5-10 keV bands respectively. 

For the EPIC imaging observing modes, photons are not only 
registered during the actual integration interval but also during the 
readout of the CCD. These out-of-time events are hence assigned 
the wrong position in the readout direction. The fraction of out- 
of-time events is highest for the PN full frame mode (6.3 %) and 
therefore for each PN exposure, an additional synthetic out-of-time 
events list was produced by randomising the coordinates of the 
events within each chip in the readout direction. Out-of-time im- 
ages were constructed in each energy band by filtering these event 
lists in exactly the same way as the real event lists. These out-of- 
time images were used as inputs to the background model as de- 
scribed in Section l2~2l 

The astrometry of the individual observations was corrected 
for small offsets between the pointings. A sourcelist was con- 
structed for each observation as described in Section l2~2l and cross- 
correlated with the optical positi ons of the 214 sources found in 
the Chandra catalogue of iMcHardv et alJ 120031) using the SAS 
task EPOSCOR. The appropriate offsets in RA and dec were then 
applied to each of the individual events to tie the XMM-Newton 
data to the optic al/Chandra /radio co-ordinate frame. The images 
and exposure maps were then reproduced with the correct astrom- 
etry. The actual offsets in RA and Dec differed slightly between 
the three observations. The first observation (revolution 276, 179 
XMM-Newton sources) had offsets of 1.4", -1.3" applied, the sec- 
ond (revolution 281, 106 XMM-Newton sources) had offsets of 
0.5", -0.5" applied and the third observation (revolution 282, 257 
XMM-Newton sources) had offsets of 0.6", 0.2"applied. 
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Rev. 


Date 




Live Time (ks) 






Filter 






MOS] 


MOS2 


PN 


MOS1 


MOS2 


276 


12.06.01 


43.1 


45.8 


35.5 


thin 


med 


281 


22.06.01 


14.1 


12.0 


6.7 


med 


thin 


282 


23.06.01 


59.2 


60.2 


47.9 


thin 


thin 



Table 1. Summary of 13 H deep field XMM -Newton observations showing 
the date and length of observations and the filters used. The live times have 
had periods of high background excluded. 



2.2 XMM-Newton Source detection 

The combined images in each energy band were source-searched 
simultaneously using the SAS tasks EBOXDETECT and EMLDE- 
TECT. EBOXDETECT is a sliding cell detection algorithm which 
outputs an initial sourcelist. This sourcelist is input for the 
EMLDETECT task which performs a maximum likelihood PSF fit 
to the sources producing refined positions and fluxes for all bands 
simultaneously. This method results in better source positions than 
searching in the individual energy bands one at a time because the 
PSF fit is based on the maximum number of counts per source. If 
the best fit source count rate in any particular energy band is less 
than zero (i.e. there are fewer counts at the source position in the 
image than in the background map) the source count rate is set to 
zero in this energy band. 

A background map was produced for each combination of ob- 
servation, detector and energy band (36 background maps in total) 
using our own software. This software performs a maximum like- 
lihood fit to the background, assuming a three-component back- 
ground model: out-of-time events, a flat unvignetted component, 
caused primarily by cosmic rays, and a vignetted component rep- 
resenting unresolved faint sources and genuinely diffuse emission. 
For the PN, the out-of-time events contribution to the background 
was fixed at 6.3% of the intensity of the synthetic out-of-time im- 
ages; for the MOS background, we assumed no contribution from 
out-of-time events. The intensities of the vignetted and unvignetted 
background components were free parameters in the fit. 

To maximise sensitivity an iterative procedure was employed. 
Initially, sources were detected in each individual image (per de- 
tector and per energy band) using a 3-pixel-square sliding cell in 
EBOXDETECT, with the background computed as the average of 
the surrounding 7x7 pixels. Then, the sources were excised from 
the individual images and the background fitted. Each of the back- 
ground maps from the MOS and PN cameras for a given energy 
band were then summed to produce one combined background map 
for each energy band. The resultant background maps were then 
used for the sliding cell (EBOXDETECT) followed by maximum 
likelihood (EMLDETECT) source detection on the combined im- 
age in each energy band. The sequence of background determina- 
tion followed by source searching was then repeated several times. 
We found that the sourcelist and background maps converged after 
4 iterations. Likelihood thresholds (DETJVIL values output from 
the source detection) of 4 and 5 were chosen for EBOXDETECT 
and EMLDETECT respectively. These values are related to the 
probability of a random Poissonian fluctuation having caused the 
observed source counts via l Cash 1979): 

DETJVIL = -lnP random (1) 

The conversion factors from count rates to flux were deter- 
mined from the EPIC response matrices, over exactly the same en- 
ergy ranges as those in which the images were constructed, assum- 



Energy (keV) 


ECF ( cts per 10 " erg cm 2 ) 


0.2-0.5 


4.7775 


0.5-2.0 


4.8905 


2.0-5.0 


1.9605 


5.0-10.0 


0.5929 



Table 2. Energy conversion factors (ECF) used to convert between count 
rate and flux. 

ing a power law spectrum with photon ind ex T _ 1.7. This is a 
good average within our flux range l P age et all2003tlMateos etlrt] 
2005); however the sources have a range of photon indices. To as- 
sess the impact of such a spread we have calculated the expected 
conversion factors using photon indices of T = 1.4 and T = 2.0 
respectively. This range represents the expected spread in spec- 
tral slope for sources contributing to the XRB. The relatively flat 
T = 1.4 lower limit corresponds to the XRB slope in the 3-15 keV 
range, produced by absorbed AGN. The upper limit value is typ- 
ically found in unabsorbed AGN, and we would therefore expect 
the bulk of our sources to lie between the two values. 

The largest effect is in the 0.2-0.5 keV and 5-10 keV energy 
bands where the conversion factors derived are different from those 
assuming a photon index of T = 1.7 by up to 11% and 8% re- 
spectively. However, in the 0.5-2 and 2-5 keV bands the photon 
index chosen only affects the conversion factors by 1-2%. The PN 
response matrices include only single and double pixel events, but 
for the three highest energy bands our PN images also include triple 
and quadruple events. In order to take this into account, the count 
rate to flux conversion facto rs were cor rected (by up to 6% in the 
hardest band) as described in Osborne et alj l2001f) . 



3 MONTE CARLO FIELD SIMULATION 

We have used an XMM-Newt on spe cific extension of the simula- 
tion method of lHasinger et alj 1199 8?) to obtain a quantitative char- 
acterisation of the source detection procedure and to assess the 
reliability of the resultant sourcelist. We have used our simula- 
tions to find the appropriate detection threshold to be applied to 
the 13" field. A Monte-Carlo approach is particularly powerful 
near the survey flux limit where a number of different processes 
contribute to uncertainties in the detected source parameters. Our 
simulation method consists of several modular steps that are re- 
peated for a large number of synthetic fields. Briefly, an 'input' 
sourcelist was generated independently in each energy band. Each 
list was then folded through the XMM-Newton imaging character- 
istics to generate images in each energy band. These images were 
then source searched to produce 'output' sourcelists. A pairing al- 
gorithm was used to associate an 'input' source with each 'output' 
source. Each of these stages are described in detail in Appendix A. 
Here we present the results of a comparison between the output and 
input source properties. These provide an indication of the biases 
inherent in our survey. One thousand simulated fields were used to 
reduce statistical uncertainties in the analysis. 

3.1 Comparison of simulated input and output sources 

In order to assess the accuracy of source positions and fluxes, and to 
estimate the degree to which confusion and Eddington bias affects 
the source counts in our 13 w data we have compared the input and 
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Figure 1. Greyscale images showing the distribution of positional offsets between output and input source locations as a function of input flux in two energy 
bands. Those sources within the central 9' of the XMM-Newton field of view having DETJV1L > 5 are shown. The concentration of positional offsets as a 
function of 5 ;„ ; , is indicated by the darkness of the greyscale image. The three contours show the distances within which 68, 90, and 95 per of the data lie. 



output properties of our simulations. The simulations should mimic 
any biases found in the real data. 

We matched each output source found in the simulated images 
to the closest valid input source. We consider an input source as 
valid when its flux, S,„ p , contributes a reasonable fraction (>20 per 
cent) to the total output flux, S BUt . No upper limit was applied to the 
radius at which input and output sources were matched in order to 
assess the typical offsets between input and output positions. Fig.Q 
shows the distribution of measured positional offsets as a function 
of input flux, Sj„ p . The greyscale image shows the relative density 
of sources at a given 5, „ p . The dark band shows where the majority 
of sources lie. All sources with DETML values > 5 and offaxis 
angles < 9'are shown. The contour lines plotted correspond to the 
positional offsets within which 68, 90 and 95 per cent of the data 
lie. 

The mean positional offset decreases with increasing flux and 
is < 10" for all but the faintest fluxes. In the 0.5-2 keV energy 
band, 95 per cent of sources with S inp > 5 x 10~ 16 erg cm" 2 s -1 and 
offaxis angles < 9' have positional offsets < 10". In the 5-10 
keV energy band, 95 per cent of sources with 5,„ p > 5 x 10~ 15 
erg cm" 2 s"" 1 and offaxis angles < 9' have positional offsets < 10". 
Those with offaxis angles > 9' have systematically larger posi- 
tional offsets (larger by ~ 2" over the majority of flux ranges). Any 
sources with higher positional offsets are most likely due to incor- 
rect associations. A discussion of the positional accuracy found in 
the real 13" data and its comparison with the simulations is de- 
ferred to Section l4~3*l 

Subsequently, we matched output sources to input sources 
within a radius r cut = 5", 8", and 10" for offaxis angles of - 9', 
9-12' and > 12' res pectively, reflecting the degradation of the 
XMM PSF tarchll 2004) away from the optical axis. Where more 
than one candidate input source lay within r cut , the brightest was 
chosen. The brightest input candidate within r a „ must be the 'cor- 
rect' input counterpart in the sense that it is the largest contributor 
to the output source counts. In practice, when averaged over the 



1000 simulations, only 0.6, 1.7, 1.5 and 0.5 output sources per field 
had more than one valid (55,„ p > S m „) input candidate within r a „ 
in the 0.2-0.5, 0.5-2.0, 2.0-5.0, and 5.0-10 keV energy bands. 

The output fluxes and positions were then compared with 
the corresponding input values. There are several reasons why we 
might expect a difference between the input and output flux distri- 
butions and source counts: i) a systematic or statistical flux mea- 
surement inaccuracy, ii) source confusion, iii) statistical fluctua- 
tions in the background which may be detected as sources, and iv) 
Eddington bias. All these factors must be considered together when 
interpreting the 13 H field data. 

The intrinsic accuracy of the source detection photometry is 
best evaluated at high fluxes and low offaxis angles, where ii), iii) 
and iv) are less important. The bright end of Fig.|5|illustrates the 
high fidelity of the detected fluxes, S out , to the input fluxes, Si„ p . 
Considering sources with 5,,,,, > 5x 10~ 14 erg cnT 2 s , with offaxis 
angles < 9' the average S oll ,/S inp ratio is 1.01±0.03 (lcr) in energy 
bands 0.2-0.5, 0.5-2 and 2-5 keV. In our highest energy band the 
average 5 D „,/5,„ P ratio is also 1.01 but the scatter is larger (lcr = 
0.06). 

The distribution of S ollt /S ,„ p is shown in Fig.[3]for three flux 



intervals in each energy band. At bright fluxes (5 



2 x 10- 14 



erg cm 2 s 1 ) the distribution is narrow, symmetrical and centred 
1, because the statistical errors on the fluxes are 



) the 



on S oul /Sj, 

small. At intermediate fluxes (S,„ p = 6 x 10~ 15 erg cirT 2 s 
distributions are still relatively symmetrical, but they are slightly 
broader as the statistical errors on the fluxes are larger. However, at 
the faintest fluxes, the distributions are much broader and signifi- 
cantly skewed towards larger S ou ,IS inp ratios. The increased width 
is due to the increased statistical errors on the fluxes. The distribu- 
tions are shifted towards larger 5 „„, / 5 mp ratios because at such faint 
fluxes sources are unlikely to be detected unless they are enhanced 
by Poisson fluctuations or by source confusion. 

Source confusion occurs when two or more nearby input 
sources fall in a single resolution element of the detector and re- 
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Figure 2. Scatter plots of input vs output flux in two selected energy bands. We only show those detected sources having DET.ML > 5 and input counterparts 
within 5", 8", and 10" for input offaxis angles of - 9', 9 - 12' and > 12' respectively. Sources with S oa ,/(S,„ p + 3<x „,) > 1.5 are plotted as triangles. The 
results for 100 simulations are shown for clarity. 
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Figure 3. The distribution in output (S out ) flux to input flux (S;„ p ) ratio for three illustrative input flux ranges. The solid line represents the sources in the flux 
interval centred on S j fl p = 2 X 10~ erg cm~ s (1 x 10 -14 - 4x 10~ 14 erg cm 2 s 1 ). The dotted line represents sources in the flux interval centred on 
Sinp — 6 x 10 erg cm - 2 s _1 (3x ur 15 - 1.2 x ur 14 erg cm s ) and the dashed line represents sources in the flux interval centred on Sj„ p = 4x10 
erg cm" 2 s" 1 (2 X 10~ 16 - 8 X lO" 16 ). The flux intervals were chosen such that AlogS,„ /; = 0.3. 
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Figure 4. a) Fraction of flux amplified sources as a function of flux, b) Fraction of unmatched sources as a function of flux. (0.2-0.5 keV (dotted), 0.5-2 keV 
(solid), 2-5 keV (dashed), 5-10 keV (dot-dashed)). The fraction of flux amplified sources is < 0.4 per cent in all energy bands and is significantly less than the 
fraction of unmatched sources. 



suit in a single output source. This results in a flux amplification in 
the output source and a net loss of fainter sources. The position of 
the output source will be close to the centroid of the merged input 
sources. Therefore when two input sources of similar flux are con- 
fused, the output position does not correspond to either of the input 
positions. Source confusion can limit the depth of any deep survey 
depending on the size of the telescope beam, and on the sky-density 
of objects as a function of flux. In practice, we cannot distinguish 
between a source boosted by photon noise or one confused with 
another faint source, therefore we consider the two effects jointly. 
We class sources as 'flux amplified' (corresponding to 'confused' 
sources in lHasinger et alll99 8) if S „,/(S,„ P + 3cr „,) > 1.5 (where 
cr oul is the ltr error on the output flux S „,„). Fig.|3t shows the frac- 
tion of flux amplified sources as a function of input flux in our four 
energy bands. The fraction is less than 0.4 per cent at all fluxes 
in all energy bands. This fraction will depend on the exact def- 
inition of flux amplification. Using a less stringent definition of 
S outliS m P + ~iCout) > 1.2 still results in a fraction well below 2% in 
each energy band at all fluxes. 

We class an output source as 'unmatched' when there are 
no valid (5 D „ f < 55 ,„„) input sources within r cut (corresponding 
to 'spurious' sources in Hasing er et aljfl993) . These are mainly 
caused by positive fluctuations in the background. Fig. 05 shows 
the fraction of unmatched sources as a function of flux in each en- 
ergy band. As expected, the unmatched fraction is highest at very 
faint fluxes, where up to 30% of the sources are unmatched. Un- 
matched sources are many times more numerous than the flux am- 
plified sources at any flux. 

In our simulations we curtailed our input sourcelists at fluxes 
5x fainter than those found in the 13 H data in each energy band 
in order to speed up processing time. In order to assess the impact 
of the simulation flux limit on the number of flux amplified and 
spurious sources, we have also produced a smaller number of sim- 
ulations to a greater depth, reaching fluxes lOx fainter than those 
found in the 13 H data. The fraction of flux amplified sources in 



these faint simulations agrees with the fraction found in our orig- 
inal simulations to within 0.02 per cent. Likewise, the fraction of 
unmatched sources agrees to within 2 per cent. We are therefore 
satisfied that our chosen flux limits are sufficiently deep. 

In order to investigate the flux limits at which confusion noise 
dominates over Poisson errors, we have produced and source- 
searched a small number of ultra-deep simulations with no Poisson 
noise. The results are presented and discussed in Appendix B, and 
show that the 13" flux limits are more than a factor of 4 brighter 
than the ultimate XMM-Newton confusion limit in any of the 4 en- 
ergy bands. 

Eddington bias lEddi ngtonll913h results in a systematic offset 
in the number of sources detected at a given flux. The magnitude 
of this effect depends on the both the statistical errors on the mea- 
sured flux values and on the intrinsic slope of the N(S). As there are 
generally many more faint sources than bright ones, uncertainties 
on the measured flux values will result in more faint sources be- 
ing up-scattered than bright sources being down-scattered. There- 
fore we would expect more faint sources to be detected than are 
actually input, and the output source counts at a given flux to be 
greater than those input. In the case of our simulations the situa- 
tion is further complicated due to the double powerlaw form of the 
N(S ) distribution and the fact that the flux error distribution is non- 
uniform and a function of several parameters including flux and 
offaxis angle. 

The level of Eddington bias expected in the 13" deep field is 
shown in Fig.|5|where the simulated input and output source counts 
are compared. Below S knee the ratio rises as the statistical errors on 
the flux measurements increase. At the lowest flux interval there 
is a drop in output source counts. The reason for this is the strong 
skew in S olll /S inp at the faintest fluxes (see Fig. [5) which boosts the 
output fluxes. The output/input source counts ratio is a minimum at 
S knee where both the statistical errors on the flux measurement are 
low and the source counts become flat. Above S k„ee the output/input 
source count ratio is constant within the errors with a value of ~ 
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Figure 5. Simulated output (crosses) 0.5-2 keV N(S) distribution nor- 
malised by the input distribution. The error bars represent Poisson errors 
on the output source counts. The dotted line represents the case where the 
input and output source counts are equal. The source counts are most dis- 
parate at faint fluxes where the output counts are enhanced by up to 23%. 



1.05 indicating that Eddington bias affects our bright counts at the 
5% level. 



3.2 Assessing the reliability of the 13" sources 

In this section we describe how we determined from our sim- 
ulations the appropriate detection threshold to be applied to the 
13 H field EMLDETECT sourcelist. Our aim was to produce a 
sourcelist such that nearly all erroneous detections are removed 
whilst retaining the maximum number of real sources. 

For each detected source, EMLDETECT measures the detec- 
tion maximum likelihood statistic, DETJV1L, which takes account 
of source counts, background counts, and the PSF. For the real data, 
the SAS task EMLDETECT provides both single-band and multi- 
band measurements of DET_ML. For the simulated images we are 
limited to single band measurements as the four bands are sim- 
ulated independently. To simulate all four bands simultaneously 
would require a priori knowledge of the sources' intrinsic X-ray 
spectra, redshifts and column densities which we do not have. The 
value of DETJVIL for any single detection is directly related to the 
probability of the source being caused by a random Poisson fluc- 
tuation via Eqn.0 However, it is difficult to translate a minimum 
threshold value of DETJVIL applied to the whole sourcelist into a 
total number of expected spurious sources in the field as the prob- 
ability is a function of position within the field, due to the varying 
PSF and exposure map. This is therefore best explored via a large 
number of Monte-Carlo simulations. 

Using the simulated sourcelists in each band, we calcu- 
lated the fraction of sources which were either flux amplified 
or unmatched as a function of the minimum detection threshold 
DETJvIL m j n . This is shown in Fig.|6| To restrict the fraction of bad 
sources in our final 13" sourcelist we use the value of DET_ML rain 
in each band at which only 5 percent of sources are either flux am- 



Figure 6. Plot showing the fraction of flux amplified or unmatched sources 
for a range of lower DET.ML limits. Results are shown for each energy 
band, 0.2-0.5keV (large dash), 0.5-2keV (solid line), 2-5keV (small dash), 
5-10keV (dot-dash). The dotted line shows the 5% badness level. 



plified or unmatched in our simulations; these are 5.9, 5.9, 6.0, 
8.1 in the 0.2-0.5, 0.5-2, 2-5, and 5-10 keV bands respectively. 
For a source to pass our significance threshold we require that it 
has DETJVIL > DETJVIL m j n in at least one energy band. In prac- 
tice, because we use a multi-band source detection process for the 
13" data, we expect fewer than 5% bad sources after we have ap- 
plied this criterion, as many of the sources will be detected in more 
than one energy band. 

For each source in the 13" sourcelist, we can identify all the 
output sources from the simulations which lie within 2' of the real 
13" source and have a similar DETJVIL (within 10%). The frac- 
tion of these output sources which are unmatched or flux amplified 
gives a good estimate of the probability that the detection of the 
real source in this energy band is unreliable. For sources detected 
in more than one band the probability that the source is spurious 
is given by the product of these individual probabilities. The total 
number of spurious sources expected in the 13" field can then be 
estimated by summing the probabilities from each source. 



3.3 Maximum likelihood N(S ) fitting method 

Our simulations are ideal for testing our N(S ) fitting procedure be- 
cause we know a priori the input N(S ) fitting parameters. Accurate 
fitting of any N(S) relation requires a knowledge of the sky area 
searched and the probability of detection at a given flux. In Fig. 
Qthe effective area of the survey as a function of limiting flux is 
shown for each energy band. This was determined by comparing 
the number of output to input sources at each flux and multiplying 
the resultant fraction by the geometric area of the survey. Addition- 
ally, the distribution of S out /Si np is required to account for Edding- 
ton bias (see Fig. [3] and Section l3~TV 

Adapting the maximum likelihood method of iMurdoch et alJ 
( 1973), we fitted a double powerlaw model to the output 
N(S) relation. Rather than fitting the observed fluxes directly, we 
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0.20 




Figure 7. Effective area for each energy band, determined from Monte- 
Carlo simulations (0.2-0.5 keV (dotted), 0.5-2 keV (solid), 2-5keV 
(dashed), 5-10 keV (dot-dashed)). 

convolved a model N(S) with the distribution of S oul /Si np to pro- 
duce a model probability distribution of observed fluxes P(S ol „). 
The best-fit N(S) was then determined by maximising the likeli- 
hood of obtaining the observed flux distribution. As a check, we 
applied this technique to our simulated data, and recover the in- 
put values to within the statistical limits of our simulations (~2%). 
We are therefore confident in applying this technique to the real 
13 H data. Fits neglecting the distribution of S out IS inp are typically 
systematically offset by ~5% 



4 RESULTS 

In this section we present our final source catalogue and results of 
fits to our source counts. We compare the properties of our sources 
with those detected by Chandra in the same area and closely exam- 
ine our quoted positional uncertainties via comparisons with both 
Chandra detections and our simulations. 

4.1 X-ray source catalogue 

Following the source detection procedure described in Section |2~21 
(simultaneous source searching in four energy bands), a total of 
275 sources were detected with DETJVIL > 5 in at least one energy 
band. This sourcelist was curtailed using the DET_ML m j n values for 
each energy band determined in Section l3~2l This reduced the final 
number of sources to 225. The final sourcelist is presented in Table 
[9] Using the procedure described in Section l3~2l we expect a total 
of 7 spurious sources. 

4.2 13" deep field source counts 

We show the integral N(> S ) in each energy band in Fig. [8] We 
have fitted single and double powerlaw models to the unbinned dif- 
ferential source counts in all energy bands using the method de- 
scribed in Section l3~3l The four identified Galactic stars in the field 



Energy (keV) y K KS Prob 



0.2-0.5 
0.5-2.0 


1 84 +a2 ° 

1,0 -0.18 

1.69+ 011 


47±20 
112+23 


1.8x10" 
1.0x10" 


-2 
-10 


2.0-5.0 




126±28 


3.5x10" 


-4 


5.0-10.0 


Z80!i 


150±18 


2.9x10" 


-5 



Table 3. Best-fit values for a single powerlaw fit to the 13 H deep 
field differential source counts. All errors are at 95% confidence. The 
best-fit slope, y, is listed together with the normalisation, K, in units of 
(10~ 14 erg cm" 2 s~ 1 ) y ~ 1 deg~ 2 . KS null-hypothesis probabilities on the fits 
are listed in the final column. 

have been excluded from this analysis. The best-fit double power- 
law models are overlaid in Fig.|8|(single powerlaw in the 5-10 keV 
energy band) as solid lines. The 95% confidence interval of the fits 
are indicated by the bowties. These incorporate errors on the slopes, 
normalisations and knee if applicable. 

The best-fit parameters for a single powerlaw fit are listed in 
Table l4.2l The uncertainties on the fit parameters are quoted at the 
95% confidence interval for one interesting parameter. We have 
tested the goodness of fit of these models using the Kolmogorov- 
Smirnov (KS) test, and the null-hypothesis probabilities that we 
obtain are given in Table|42l The fits are unacceptable in all bands. 

For the double powerlaw fits we have fixed the slope at bright 
fluxes (72) and fit only for the slope at faint fluxes (71) and knee 
position (S k,we)- In the 0.2-0.5 keV band we set 72 = 2.51 based 
on results from the ROSAT All Sky Survey IVoees et aljl999t) . The 
value of 72 in the 0. 5-2 keV band was fixed at 2.60, derived from 
the RIXOS survey iMason et al1l2000l) . For the 2-5 keV band 72 
was fixed at 2.65 as found in the BeppoSAX High Energy Large 
Area Survey (HELLAS iGiommi et alJl200d). and consistent with 
the ASCA derived value of 2.5 from lUeda et alj ll999l) . The fit pa- 
rameters for the double powerlaw fits are listed in Tabled The 
double powerlaw model provides a better fit compared to a single 
powerlaw model in each of these three energy bands. However, the 
double powerlaw model fit is formally rejected with 99% confi- 
dence in the 0.5-2 and 2-5 keV energy bands. In all three energy 
bands S k„ ee occurs between 1.08 - 1.27 x 10 -14 erg cm" 2 s" 1 . 

In the 5-10 keV band the slope found in the singl e powerlaw 
fit is consistent with that found by iBaldi et alJ (2002) at brighter 
fluxes. We attempted to fit a double powerlaw model to the 5-10 
keV band source counts, but in this case 5fa, ee is unconstrained. 
We therefore consider that a double powerlaw model for the source 
counts is not justified at the depth of our survey in this energy band. 

4.3 Positional uncertainties 

To assess the reliability of the EMLDETECT positional errors in 
the XMM-Newton sourcelist we have cross-correlated the final 
XMM-Newton and Chandra sourcelists. The spread in positional 
offsets between the two catalogues should provide a good repre- 
sentation of the true spread in XMM-Newton positional errors. To 
prevent the (albeit small) statistical errors on the Chandra posi- 
tions contributing to the spread, where possible we used the po- 
sition of the optical counterpart to the Chandra source: the po- 
sitiona l errors on the optical counterparts are typically less than 
0.3" JrvIcHardv et ai]|2003l) . The two sourcelists were matched 
within a search radius of 10" over the entire XMM-Newton field 
of view. 

The two catalogues have 155 sources in common. The posi- 
tional offset between the 155 matched XMM-Newton and Chandra 
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Figure 8. 13 deep field integr al N(> 5) in each energy band. Bowties indicate 95% errors. O verlaid for compa rison are the results from the 
CDF-S (dashl iRosati et all l2002h CDF -N (long dash) i Bran dt et alJl200ll) . HELLAS2XMM (dot) iBaldi et alJl2002h and ASCA f dot-dot-dot-dashl 
iCagnoni, Delia Ceca & Maccacaro 1998) surveys. In the 0.5-2.0 keV band the ROSAT counts found in the 13" field <McHardv et alll998l) are also overlaid 
(dot-dash). In the 2-5 keV band the results from Cowie et al. l2002tl are also overlaid (dot-dash). Triangles in the 5-10 keV band denote the BeppoSAX counts 
o flFiore et ai]|200ll) . 
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Table 4. Best-fit parameters for a double powerlaw fit to the 13 H differential source counts. The slopes and normalisations below and above the break flux, 
Sknee are denoted by i and 2 respectively. The values of yi in the three energy bands were fixed to appropriate values found from the literature. The break flux, 
Skme is in units of 10~ 14 erg cnT 2 s~'. Normalisations K\ and K2 are in units of (10~ 14 erg cnT 2 s~ [ ) y ~ l deg~ 2 . All errors are at 95%. KS null-hypothesis 
probabilities on the fits are listed in the final column. 
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Figure 9. Positional offsets between XMM-Newton and Chandra counter- 
parts as a function of XMM-Newton offaxis angle. Triangles denote XMM- 
Newton detections with > 1 Chandra counterparts within 10" (numbers 
indicate the XMM-Newton flux in units of 10~ 14 erg cm -2 s - ' ). The shaded 
areas indicate the regions where 68% (dark grey) and 95% (light grey) of 
the simulated data lie. In the case of the simulations the positional offset is 
that between the input and output source. 86% of the 13" data lies within 
the 68% simulation contour, indicating that the real data have better source 
positions. This is due to the fact that the real data are generally detected in 
> 1 energy band. 



sources as a function of XMM-Newton offaxis angle is shown in 
Fig. [5] The positional offsets increase with XMM-Newton offaxis 
angle, but 80 per cent of the matched sources have positional off- 
sets < 2"; 95 per cent of the sources have positional offsets < 4". 

We compared this distribution with that found from the sim- 
ulations (see Section ITTl which should provide a good indication 
of the distribution of positional offsets arising within the source de- 
tection chain. The contours in Fig.[9]show the distribution of differ- 
ences between input and output source positions in the 0.5-2 keV 
simulations. The dark (light) grey area shows where 68% (95%) 
of the simulated sources lie. The distribution of XMM-Newton - 
Chandra offsets is more highly peaked than the input-output off- 
sets in our simulations: 86% of XMM-Newton -Chandra offsets lie 
within the 68% simulation contour. This is because the real sources 
in the 13" field are generally detected in more than one energy 
band, contrary to the simulated sources. Fitting the source positions 
in all four energy bands simultaneously, as was done to produce our 
source catalogue, results in better positional accuracy than a single 
band fit because more source counts are used in the fit. 

Five of the matched sources have more than one Chandra 
counterpart. However, in two out of the five cases, the closest coun- 
terpart is significantly brighter than the other one, and we are easily 
able to identify the correct counterpart. This suggests that ~ 2% of 
our sources are confused, in broad agreement with the fraction of 
flux amplified sources expected from our simulations. 

In Fig.[l0| we plot the statistical positional uncertainties of the 
XMM-Newton sources given by EMLDETECT against the actual 
positional uncertainties given by the XMM-Newton - Chandra off- 
sets. These are compared with our 0.5-2 keV simulation distribu- 
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Figure 10. Positional offsets between XMM-Newton and Chandra coun- 
terparts compared with quoted XMM-Newton positional error. The shaded 
areas indicate the regions where 68% (dark grey) and 95% (light grey) of 
the simulated data lie. In the case of the simulations the positional offset is 
that between the input and output source. 43% of the 13"data lie within the 
68% simulation contour indicating that EMLDETECT underestimates the 
true positional error. 

tions where in the case of the simulations the positional offsets are 
those of the input/output source positions. The dark (light) grey ar- 
eas show where 68% (95%) of the simulated sources lie. In general, 
the positional offsets are slightly larger than the EMLDETECT er- 
rors for the 13" data. However, for EMLDETECT positional errors 
less than 1" the positional offsets have a broader distribution in the 
13" data than in the simulations: only 43% of the 13" data lies 
within the 68% simulation region. This implies that the EMLDE- 
TECT positional error in the 13" field underestimates the true 
positional error. We have therefore added in quadrature a system- 
atic error to the statistical EMLDETECT positional error of each 
13" source, such that 68% of the 13" sources lie within the 68% 
region of the simulations. This is achieved with an additional sys- 
tematic positional error of 0.76" (l<x). This additional component 
may be due to residual uncertainties in the detector geometry and 
may represent a fundamental limit to the accuracy of any XMM- 
Newton position. According to our simulations, EMLDETECT po- 
sitional errors greater than ~ 2.5" may underestimate the true po- 
sitional error by an order of magnitude (the 95% simulation region 
limit samples a non-gaussian tail at this point) and should be used 
with this caveat in mind. Section l5~2l further discusses the reliability 
of the survey and the implications of the positional uncertainties. 

4.4 Comparisons with Chandra detections 

There are 70 XMM-Newton sources with no Chandra counterpart. 
Of these, 9 sources are either on t he very edge or out of the field of 
view of the Chandra mosaic (see lMcHardv et alj2003l) and would 
therefore not be expected to have a Chandra counterpart. R-band 
images, centred on the positions of these 70 sources are presented 
in Fig.EI] (Optical im ages of those so urces with a Chandra coun- 
terpart are presented in McH ardvetaTT2 003). In 49 cases there is a 
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Figure 11. Top: 10" X 10" i?-band Subaru SuprimeCam images centred on the XMM-Newton sources without Chandra counterparts. The XMM-Newton 
reference number (shown in column 1 of Table l9l is at the top left of each image. The XMM-Newton off-axis angle is at the top right. The XMM-Newton 
0.2 - 10 keV flux (in units of 10~ 15 erg cirT 2 s _1 ) is at the bottom left and the optical counterpart R magnitude, if applicable, is given in the bottom right. 
All of these Suprim eCam imag es have the same greyscale levels to ease comparison of the optical counterparts' brightnesses. The greyscale was chosen to 
correspond to that in McHardv et al. 12003). Bottom: 30" X 30" i?-band Subaru SuprimeCam images centred on the 6 XMM-Newton sources with Chandra 
counterparts which are too extended to fit in a 10" X 10" box. Source number 18 is actually shown in a 70" X 70" box. The top images are at the same grey 
scale level as Fig. l 111 The lower images have their grey scales adjusted for maximum contrast. Source number 10 is not listed as having an optical counterpart 
in Tablel5las the centre of the barred spiral galaxy is outside the error circle of the XMM-fctoii source. However, we may be seeing X-ray emission from an 
X-ray source in the spiral arms of the galaxy. 
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XMM-Newlou 


RA J(2000) 


Dec J(2000) 


Pos. err. 


Off axis 


0.2-10.0 keV 


Extension 


R mag. 


Chandra 


Faint Chandra 


Comments 


No. 






(arcsec) 


(arcmin) 


Flux 






counts 


source? 




1 


13 33 28.96 


37 55 58.82 


1.57 


12.91 


14.86 


- 


19.74 


3 


n 




2 


13 33 30.32 


37 56 12.21 


2.38 


12.67 


5.94 


- 


23.18 


1 


n 




3 


13 33 30.68 


37 55 05.57 


3.04 


12.53 


2.55 


- 


25.83 





n 




10 


13 33 38.91 


38 01 56.15 


4.83 


12.96 


7.88 


- 


0.00 





n 




1 ] 


13 33 40.70 


37 49 47.79 


1.99 


11.74 


1.65 


- 


23.76 


2 


n 




18 


13 33 43.80 


37 54 55.76 


2.56 


9.95 


1.74 


- 


15.63 


3 


y 




20 


13 33 45.45 


37 58 08.04 


0.92 


10.14 


8.82 


- 


17.45 





n 




22 


13 33 47.52 


37 53 51.51 


1.11 


9.27 


3.46 


- 


23.39 


1 


n 




26 


13 33 50.80 


37 57 16.61 


1.77 


8.88 


4.43 


- 


25.41 


3 


y 




28 


13 33 54.21 


38 02 53.17 


1.23 


11.21 


16.49 


- 


25.41 





n 




31 


13 33 57.90 


37 49 59.01 


2.91 


8.70 


0.98 


- 


24.79 


1 


n 




33 


13 33 58.87 


38 00 24.88 


1.18 


8.87 


6.60 


- 


24.83 


1 


n 




35 


13 33 59.56 


37 51 39.75 


2.34 


7.57 


2.06 


- 


25.69 


2 


n 




36 


13 33 59.63 


37 49 26.35 


12.19 


8.75 


1.86 


- 


17.64 





n 




40 


13 34 01.06 


38 00 25.75 


1.50 


8.55 


5.92 


- 


20.50 


1 


n 




44 


13 34 07.18 


37 59 59.34 


1.23 


7.36 


4.87 


- 


22.50 





n 




47 


13 34 08.39 


37 47 50.85 


2.69 


8.72 


4.67 


- 


22.65 


1 


n 




54 


13 34 12.36 


37 59 10.36 


2.01 


6.06 


1.10 


- 


26.57 


1 


n 




56 


13 34 13.59 


37 45 39.06 


1.22 


10.12 


2.42 


- 


16.43 


3 


y 




62 


13 34 15.49 


38 03 05.78 


1.80 


8.97 


4.81 


- 


19.68 


2 


n 




63 


13 34 15.54 


37 52 27.80 


1.31 


4.43 


5.54 


- 


23.65 


2 


y 




66 


13 34 17.80 


37 44 31.39 


2.04 


10.89 


3.51 


- 


18.76 


2 


y 




70 


13 34 19.56 


37 51 47.77 


1.40 


4.26 


3.58 


- 


22.89 





y 


Possible faint source outside error circle 


79 


13 34 22.70 


37 55 23.49 


2.15 


2.32 


0.9(1 


- 


18.17 


3 


n 




82 


13 34 25.91 


37 54 59.90 


1.58 


1.64 


1.63 


- 


23.94 


2 


n 




83 


13 34 26.87 


38 00 27.92 


1.45 


5.73 


2.66 


- 


0.00 


1 


n 




87 


13 34 28.66 


37 57 48.26 


1.92 


3.09 


2.21 


- 


20.64 


3 


y 




88 


13 34 28.80 


37 53 37.96 


1.78 


1.67 


2.15 


- 


18.46 





n 




89 


13 34 29.20 


38 02 45.33 


1.79 


7.90 


0.S9 


- 


22.56 





n 




90 


13 34 29.30 


38 06 51.45 


1.76 


11.98 


7.29 


- 


21.88 


1 


n 


Edge of Chandra FOV 


1(14 


13 34 34.92 


38 07 03.03 


1.99 


12.13 


6.71 


- 


24.47 


3 


n 




1 10 


13 34 36.38 


38 05 14.02 


1.25 


10.32 


5.16 


- 


25.32 


4 


y 




1 16 


13 34 38.83 


37 40 22.03 


1.11 


14.58 


34.78 


- 


25.30 





n 


Outside Chandra FOV 


117 


13 34 38.96 


37 42 43.95 


1 .60 


12.22 


5.49 


- 


24.30 





n 




118 


13 34 39.75 


37 57 01.93 


1.29 


2.38 


2.91 


- 


16.14 


4 


y 




120 


13 34 42.29 


37 41 46.46 


0.78 


13.24 


10.84 


- 


20.14 





n 


Outside Chandra FOV 


122 


13 34 42.81 


37 42 42.62 


1.55 


12.33 


9.25 


- 


0.00 





n 


Edge of Chandra FOV 


123 


13 34 42.84 


38 03 53.27 


2.33 


9.13 


1 .65 


- 


18.53 


1 


n 




133 


13 3447.18 


37 57 14.88 


1.53 


3.46 


3.71 


- 


23.73 


4 


y 




136 


13 34 48.22 


37 54 14.33 


1 .06 


2.85 


7.93 


- 


0.00 





n 




139 


13 34 51.27 


37 40 51.30 


1.23 


14.46 


17.14 


- 


0.00 





n 


Outside Chandra FOV 


141 


13 34 51.86 


38 06 35.33 


1.78 


12.18 


4.16 


- 


22.80 


1 


n 




144 


13 34 53.13 


37 48 18.85 


1.76 


7.59 


3.59 


- 


24.06 


2 


y 




147 


13 34 54.55 


37 42 07.32 


1.30 


13.41 


12.44 


- 


0.00 


4 


n 




149 


13 34 55.11 


37 49 51.41 


2.44 


6.53 


2.87 


- 


25.26 


1 


n 




150 


13 34 55.41 


37 50 40.54 


2.66 


5.96 


1.32 


- 


24.95 


4 


n 




151 


13 34 56.39 


37 39 54.75 


2.06 


15.63 


15.15 


- 


24.89 





n 


Outside Chandra FOV 


153 


13 34 56.71 


37 50 25.09 


1.74 


6.32 


4.97 


- 


0.00 





n 




154 


13 34 56.76 


37 56 03.31 


1.48 


4.59 


3.04 


- 


20.66 





n 




155 


13 34 56.87 


37 52 48.73 


2.34 


4.94 


0.63 


- 


25.00 





n 




159 


13 34 58.11 


37 57 33.15 


1.22 


5.40 


3.34 


- 


0.00 





n 




160 


13 34 58.14 


37 50 52.22 


2.32 


6.22 


1.16 


- 


0.00 


1 


n 




163 


13 34 59.48 


37 57 39.65 


1.17 


5.69 


7.07 


- 


26.38 


2 


n 




172 


13 35 05.13 


37 45 27.22 


1 .90 


11.26 


3.52 


- 


0.00 





n 




177 


13 35 07.69 


37 48 27.88 


2.66 


9.24 


2.74 


- 


25.23 





n 




179 


13 35 09.25 


38 04 01.91 


1.55 


11.43 


6.08 


- 


0.00 


1 


n 




182 


13 35 12.27 


37 48 54.89 


1.86 


9.62 


1.73 


- 


24.32 





n 




185 


13 35 14.41 


37 49 12.11 


1.49 


9.78 


20.46 


9.2±0.4 


20.55 





n 




194 


13 35 17.03 


37 49 13.98 


14.06 


10.19 


1.44 


- 


19.44 


15 


n 


Large error circle 


2(15 


13 35 24.72 


37 51 24.06 


1.68 


10.57 


3.14 


- 


0.00 





n 




206 


13 35 24.88 


37 44 38.51 


1 .90 


14.34 


3.94 


- 


23.41 





n 




209 


13 35 32.22 


38 03 39.97 


0.70 


14.39 


238.21 


62±104 


25.10 





n 




210 


13 35 32.74 


37 45 14.16 


4.66 


15.08 


71.49 


26+2 


22.43 


2 


n 




21 1 


13 35 33.05 


37 48 02.14 


1.98 


13.50 


3.24 




18.54 


1 


n 




215 


13 35 37.26 


37 47 22.80 


1.34 


14.55 


5.96 




20.28 





n 




217 


13 35 39.37 


37 56 14.96 


0.68 


12.92 


20.15 




20.89 





n 


Outside Chandra FOV 


218 


13 35 42.49 


37 55 42.93 


0.59 


13.49 


23.77 




19.47 





n 


Outside Chandra FOV 


219 


13 35 43.17 


37 52 58.65 


1.72 


13.74 


6.49 




26.58 


1 


n 


Chandra 102 is ~ 10" North 


220 


13 35 43.26 


37 53 24.27 


2.17 


13.70 


3.95 




25.89 





n 


Chandra 102 is - 15" South 


225 


13 35 58.70 


37 54 06.31 


2.76 


16.68 


17.44 


6.8±1.2 


17.51 





n 


Outside Chandra FOV 



Table 5. Properties of XMM-Newton sources not found in the Chandra catalogue of iMcHardv et all l20O3t) . The XMM-Newton flux is in units of 10~ 15 
erg citT 2 s , The extension of the source (if applicable) in the XMM-Newton image in units of arcsec. The column marked 'Chandra counts' gives the 
number of counts detected within the XMM-Newton error circle in the Chandra image. Excess emission associated with a possible faint source is indicated by 
a 'y' in column 10. See the main text for further details. 
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Chandra 


XMM-Newton offaxis angle 


0.5-7 keV 


Faint XMM-Newton 


Comments 


No. 






Flux 


source? 




82 


7.0 




1.16 


y 


y away rrom Chandra 28/ XMM-Newton 50. 


91 


12.5 




0.94 


y 


Below final threshold. 


96 


13.0 




0.91 


y 


18" away from Chandra 11 {XMM-Newton 99. 


102 


13.0 




0.80 


y 


Between XMM-Newton 219 and XMM-A'ewfo« 220. 
(ll"and 15"away respectively). 


1 15 


11.8 




0.70 


y 


7"away from Chandra 107/XMM-Newton 214 


121 


10.4 




0.66 


n 




145 


10.3 




0.53 


y 


9 away trom Chandra oilXMM-Newton I7l. 


149 


7.1 




0.50 


y 


I5"away from Chandra 23PiMM-Newton 39. 


154 


11.7 




0.47 


n 




158 


11.5 




0.46 


y 


13"away from Chandra \ \7 /XMM-Newton 78. 


159 


6.5 




0.45 


y 


ll"away from Chandra 20PiMM-Newton 175. 


163 


12.1 




0.43 


y 




164 


10.9 




0.43 


n 




171 


8.6 




0.37 


n 




181 


3.9 




0.33 


n 




182 


11.6 




0.32 


y 


Below final threshold. 


183 


7.4 




0.31 


y 


T1 1 J2 1 . 1 111 

Below final threshold. 


185 


7.0 




0.31 


y 




188 


7.0 




0.29 


y 


14" away from Chandra 52/XMM-Newton 43. 


192 


4.3 




0.27 


y 


Below final threshold. 


193 


2.4 




0.27 


y 


Below final threshold. 


195 


8.0 




0.27 


n 




197 


7.5 




0.26 


y 


T-» 1 J2 1 . 1 111 

Below final threshold. 


198 


5.5 




0.26 


n 




201 


3.9 




0.24 


n 




203 


3.8 




0.23 


y 


12"away from Chandra 108/XMM-A'evvfofi 103. 
In a region of extended emission. 


204 


7.4 




0.23 


y 


18"away from Chandra WfXMM-Newton 119. 


205 


3.3 




0.23 


n 




206 


7.7 




0.23 


y 


7"away from Chandra \42PiMM-Newton 128. 


207 


7.8 




0.22 


n 




208 


7.4 




0.21 


n 




210 


1.0 




0.19 


y 


Below final threshold. 


211 


9.5 




0.19 


y 




213 


7.2 




0.19 


y 


40" away from Chandra 28/XMM-Afevv/on 50, 



Table 6. List of Chandra sources with no XMM-Newton counterpart. The Chandra flux is in units of 10~' 5 erg cirT 2 s -1 . Excess emission associated with a 
possible faint source is indicated by a 'y' in column 4. The comment 'Below final threshold' indicates a source which was detected with EMLDETECT with 
DET_ML> 5 but subsequently excluded from the final sourcelist according to the criteria in Section l3*2l 



clear optical counterpart or counterparts within the positional error 
radius of the source. Those sources which have optical counter- 
parts extended beyond 10" are shown in wider images covering 
30" x 30" below the main image. Those in the top panel have the 
same greyscale applied as in the main image, however in the lower 
panel the greyscale has been individually adjusted for greater clar- 
ity. Searching out to a radius of 3" the number of sources with 
possible optical counterparts is increased to 58. Of the 9 sources on 
the edge of, or outside, the Chandra field of view, 7 have an optical 
counterpart. 

Table [5] lists the basic properties of the 70 XMM-Newton 
sources without Chandra counterparts. The Chandra mosaic im- 
age was visually inspected at the positions of each XMM-Newton 
source to see if any faint X-ray emission could be seen that was 
not formally detected in the Chandra source searching. Details of 
the number of counts observed within the XMM-Newton positional 
error circle for each source are listed in column 10 and further com- 
ments are listed in column 11. Of the 61 sources not formally de- 
tected within the Chandra field of view, faint emission is visible in 



1 1 cases. None of the extended XMM-JVwton sources are detected 
with Chandra. The XMM-A'ew;on observations are more sensitive 
to faint extended sources than the Chandra observations. One ex- 
ample is the ROSAT source R 1 1 7 ( McHard v~t all 19981) which was 
detected with XMM-Newton (number 56) but missed by Chandra 
. This source is a faint starburst gala xy with ext ended emission, 
probably on the scale of the galaxy jGunn et all200ll) . 

Conversely, there are 34 Chandra sources with no XMM- 
Newton counterparts (within the XMM-Newton FOV). The prop- 
erties of these sources are summarised in Table |S| Seven of 
these sources were originally detected in our initial XMM-Newton 
sourcelist, but removed from our final sourcelist after we applied 
the DETJVIL cutoffs derived in Section l3~2l Visual inspection of 
the XMM-Newton images at the remaining 27 Chandra source po- 
sitions, suggests a faint source in 16 cases. As well as the 4 close 
(< 10") pairs of sources found with Chandra there are also 9 faint 
sources which were missed with XMM-Newton because they are 
in the wings of much brighter sources. This highlights the impor- 
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tance of our Chandra coverage in order to obtain accurate source 
positions. 



5 DISCUSSION 

5.1 Source density in the 13" field 

Our 13 w field represents one of the deepest blank field surveys with 
XMM-Newton and as such it is ideal for the investigation of faint 
X-ray source counts. Further, this work is the only XMM-Newton 
based survey published to date to incorporate detailed Monte-Carlo 
modelling of the detection process allowing accurate fitting of the 
source counts. Our source counts are inconsistent with both a single 
and double powerlaw fit in all energy bands, except for our lowest 
energy band (0.2-0.5 keV) where a double powerlaw fit is accept- 
able. However, a double powerlaw fit provides a better represen- 
tation of the source counts according to the KS test in all but our 
hardest (5-10 keV) energy band where we are unable to constrain 
the break. 

Previous studies show a consensus in shape: the measured 
source counts are described by a double powerlaw distribution flat- 
tening below a flux of ~ 1 x 10~ 14 erg crrT 2 s and steeper than 
an Euclidean slope at fluxes above the break (5 > S kn ee)- However 
the measured normalisat i ons differ by up to ~ 30% iBrandt et alJ 
l200ltlRosati et alJl2002T) . lYang et alj 120031) attribute this field-to- 
field variation to the fact that the fields studied cover small areas 
and are subject to cosmic variance. Here the underlying clustered 
large scale structure imprints a variation in source density on small 
scales. 

The slope of our soft band (0.5-2 keV) source counts, be- 
low the break flux of 1.08j,™ X 10~ 14 erg c m -2 s -1 , is c onsis- 
tent with recent determinations using Chandra iRosati et all2002t 
iBauer et alj|2004 lYang et alj|2004 . Our 13" field is the rich- 
est X-ray blank field observation reported to date in this energy 
band. Our overall normalisation is higher, though consistent with, 
the CDF- N results oflBrandt et alJfeOOll) and the HELLAS 2XMM 
counts of Baldi et al. (2002). However, the CDF-S source counts of 
IRosati et al] |20 02|)rand the observed counts in the Lockman hole 
( Hasi nger et alll998Ll200lh are lower by -3 0%. Allowing for th e 
errors on the CDF-S source counts, given by Rosati et al ] <2002l) . 
and Poisson errors on our counts, our field is still inconsistent with 
the CDF-S at the 2.7cr level. It is likely therefore that cosmic vari- 
ance has an important effect on the measured source counts found 
in deep X-ray fields. 

In the 2-5 keV energy band our observed counts are consi stent 
with previous Chandra and XMM-Newton stu dies jBrandtetajJ 
| 200lt lHasinger et alJl20oH IRosati et alJl2002t iBaldi et alj|2002t 
Icowie et all200llYang et all2004 . Our brightest source is a fac- 
tor of ~ lOx fainter than the limiting fl uxes of the se rendip- 
itous ASCA surv eys jCagnoni. Deli a Ceca & Maccacard Il998t 
lUeda et alj[l999l a) which makes comparisons difficult. However, 
the ASCA fits of lCagnoni. Delia Ceca & Maccacarcl JT998) appear 
to be consistent with ours when extrapolated to the flux range cov- 
ered by our survey. Our overall normalisation at faint fluxes is 
closer to the CDF-S normalisation than it is to that found in the 
CDF-N in this energy band. However, the CDF-S counts drop far 
more rapidly than ours towards brighter fluxes. 

In our hardest energy band (5-10 keV) the observed 
counts are consistent with a Euclidean slope; y = 2.80^S«. 
The slope agrees within the e rrors w ith the re s ults o f 
ICagnoni. Delia Ceca & Maccacard <1998h . iFiore et alj <200ll) . 




0.1 1.0 
S (10"" erg cm" 2 s M ) 

Figure 12. 0.5-2 keV energy band differential N(S ) distribution. Overlaid 
is our double powerlaw fit (solid line). Significant excesses over the model 
distribution occur at a flux of ~ 2 X 10~ 15 erg cm~ 2 s -1 . Also shown is the 
best fit model fro m the Chandra Survey of the Lockman Hole-Northwest 
lYang et all2004l) . which covers a similar flux range as our survey (dotted 
line). 

IBaldi et alJ 120021) and the XMM-Newt on 5-10 keV counts in 
the Lockman hole l Hasi nger et"ail 1200 lh . We expect a break in 
the source counts at a flux of ~ 4 x 10~ 15 erg cn T 2 s" 1 as re- 
porte d from Chandra observations of the CDF-S iRosati et alJ 
2002). It is therefore unsurprising that it is not detected in our 
survey, as we have no sources below this flux in our 5-10 
keV sourcelist. At the brightest fluxes our counts are in agree- 
ment with findings from the HELLAS2XMM survey iBaldi et all 
l2002h. fieppoSAX counts iFiore et alj|200ll) and ASCA counts 
ICagnoni. Delia Ceca & Maccacarcll99i}) l 

Although improved upon a single powerlaw, the double pow- 
erlaw fit to the source counts is formally rejected in all but the soft- 
est energy band (0.2-0.5 keV). The differential counts show more 
clearly where the deviations from the model occur and are illus- 
trated in Fig. ^| and Fig. 1131 There is an excess over the fit at a 
flux of ~ 2 x 10~ 15 erg cnT 2 s _1 in both the 0.5-2 keV and 2-5 keV 
energy bands. The fact that this excess is seen in both energy bands 
suggests that the feature is most likely due to clustering in the field 
rather than the X-ray spectral properties of the sources around this 
flux. We defer a study of clustering to a later paper lLoaring et alJ 
2005). 

We have also compared our source counts with wider area 
Chandra surveys which we expect to provide a more accurate value 
of the global source counts. In the 0.5-2 keV energy band our fits to 
the differential counts are simila r to the fits o btained for the Lock- 
man Hole-Northwest field (Yangetal. 2004). This implies there- 
fore that the CDF-S is a particularly underdense region. In the 2- 
5 keV energy ba nd our double powerlaw fit lies above the fit of 
ICowie et alJ i2002T) which was obtained from the source counts in 
the CDF-S, CDF-N, and two Hawaii survey fields (SSA13 and 
SSA22) observed with Chandra. If the CDF-S data is excluded 
from their analysis their faint end normalisation agrees very well 
with ours, again suggesting that the CDF-S is underdense. This 
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0.1 1.0 

S (10"" erg cm" 2 s~ 1 ) 

Figure 13. 2-5 keV energy band differential N(S) distribution. Our best 
fit double powerlaw fit is overlaid (solid line). Significant excesses over 
the model distribution occur at a flu x of ~ 3 X 10~ 15 a nd ~ 2 X 10 -14 
erg cirT 2 s _I . The best fit model of ICowie et al] 120021) is also overlaid 
(dotted line) which is based upon fits to the combined Chandra deep field 
counts and Hawaii survey fields SSA22 and SSA13. 

appears to contradict the earlier comparisons between the integral 
2-5 keV band counts in wh ich our fain t end normalis ation closely 
matches the fits from Rosati et al. ( 2002). However, the Rosati et al. 
( 2002) CDF-S counts clearly drop below ours at fluxes > 5 x 10" 15 
erg cirr 2 s -1 and it i s at fluxes above th e knee that our counts are 
most discrepant from lCowie et alJl2002l) . The comparisons depend 
upon the adopted conversion between instrument energy bands and 
also the range over which fitting is performed. Rosati et al. 1 2002) 
fit to much fainter fluxes than ICowie et alJ l2002t) . As the source 
counts flatten at faint fluxes this would result in an apparently 
higher normalisatio n at brighter fluxes and may exp lain why the 
iRosati et alJ j2002h normalisation is higher than the ICowie et alJ 
(2002) normalisation at fluxes around 2 x 10~ 15 erg ctrr 2 s _I . 

In converting values from the literature into our energy bands 
we have used a photon index of T = 1.7 which represents the av- 
erage slope of the XRB in the 0.2-12 keV energy range. Other 
authors have chosen to eit her use indi vidual photon indices mea- 
sured fr om their sourc es ( Bran dt et alj|2001l) or a relatively flat 
T = 1.4 iRosati et all2002l) which represents the overall XRB spec- 
tral slope. lCowie et alJ 120021) use a photon index of T = 1 .2. as they 
assume that the absorbed population must be significantly harder 
than the unabsorbed population to produce the overall XRB spec- 
trum. However, this photon index appears to be inappropriate at the 
fluxes probed by our survey. Consider ing the results from our X-ray 
spectroscopy iPaae et all2003ll2005l) . the vast majority of sources 
have X-ray spectra with a softer photon index than T = 1.2, al- 
though they may be absorbed in the soft band. In fact the value 
of the photon index used has little effect on source counts. Using 
a photon index of 1.2 and 1. 7 respectively to convert the 2-8 keV 
counts of Cowie et al. l2002h into our 2-5 keV energy band results 
in a integral normalisation difference of ~ 13% at our faintest fluxes 
(2.2 x 10~ 15 erg cirT 2 s _1 ): counts derived assuming T = 1.7 are 
higher than those assuming T = 1.2. This direct conversion does 



R magnitude 


Expected number 


Observed number 


18-19 


0.5 


10 


19-20 


1.1 


11 


20-21 


2.4 


18 


21-22 


5.3 


14 


22-23 


11.7 


61 


23-24 


26.1 


54 


24-25 


57.7 


42 


25-26 


127.7 


34 



Table 7. The number of XMM-Newton optical counterparts as a function 
of R magnitude. X-ray - optical cross matching was carried out within the 
95% positional error circle of each XMM-Newton source and the number of 
optical detections are displayed as a function of R magnitude. Also shown 
is the number of expected chance coincidences within e ach R magnitude 
interval given the optical source counts of lMcHardv et alJ 120031) . 



not account for the fact that the original measured fluxes would 
have been derived from count rates using different photon indices. 
We have investigated this effect further using spectra which were 
simulated in XSPEC using the Chandra ACIS-I response matrix. 
Assuming photon indices of T = 1.7 and T = 1.2 respectively to 
convert count rates to fluxes, we obtain a flux ratio F1.7/F1.2 of 0.9 
in the 2-5 keV energy band. Incorporating this factor into our flux 
conversions reduces the apparent normalisation difference to ~3% 
at a flux of 2.2 x 10~ 15 erg cirT 2 s (with the normalisation of 
the T = 1.7 source counts higher than that of the T = 1.2 source 
counts). This factor is not large enough to account for the mea- 
sured differences in normalisation in the 2-5 keV energy band. It 
is not necessary to employ any flux conversions for the 0.5-2 keV 
energy band. However, there is still a ~ 30% discrepency in source 
counts between the deep fields surveyed to date. One concern is the 
cross-calibration between the different instruments used. However, 
ASCA and Chandra fluxes agree at about 10% and Chandra and 
XMM-Newton fluxes agree at the 5% level, and so the differences in 
normalisation cannot be attributed to instrument calibration offsets 
(Snowden 2001). In the 0.5-2 keV band, we are confident that the 
field-to-field variations observed are real and due to cosmic vari- 
ance. 

We find a larger field-to-field variation in the soft band (0.5-2 
keV) source counts than in the hard band counts (2-5 keV). These 
findin gs are at o dds with the recent clustering measurements of 
lYang et al] 120031) who find hard band sou rces to be more st rongly 
clustered than soft band sources. Further, iGilli et alJ <2005h have 
found no significant difference in clustering strength between soft 
and hard sources. At present the evidence for any variation in clus- 
tering strength due to X-ray hardness appears inconclusive. Clearly 
larger datasets are required to investigate this issue further. Com- 
parisons with wide area XMM-Newton and Chandra serendipitous 
surveys coupled with further blank field observations should pro- 
vide sufficient signal to make significant advances in this area. 



5.2 Survey reliability and capabilities 

Our survey was specifically designed to study the X-ray source 
population within one decade in flux either side of the break in the 
source counts. Contrary to recent deep Chandra surveys we aim 
not to resolve the entire XRB but rather to determine the domi- 
nant emission mechanisms and amount of absorption in these faint 
sources. To this end we require high quality X-ray spectra, in par- 
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Radio Flux (mJy) 


Expected number 


Observed number 


50-100 


3.1 X 10 


1 


10-50 


4.4 x 10~ 2 





5-10 


1.8 x ur 2 





1-5 


0.3 


3 


0.5-1 


0.1 


4 


0.1-0.5 


3.3 


18 


0.05-0.10 


2.0 


6 



Table 8. The number of XMM-Newton radio counterparts as a function 
of radio flux. Radio - optical cross matching was carried out within the 
95% positional error circle of each XMM-Newton source and the number 
of radio detections are displayed as a function of radio flux. Also shown is 
the number of expected chance coincidences given the best-fit radio source 
count model of Sey mour et all 120041) . 



ticular at hard energies to probe the absorbed population. The high 
throughput of XMM-Newton at hard energies makes it particularly 
suited to this kind of study. 

To efficiently study the faint source population we require a 
highly reliable source catalogue with minimal spurious source con- 
tamination. We have examined the quality of our final catalogue 
via extensive Monte-Carlo simulations (Section[3j. Employing the 
likelihood cutoffs in each energy band determined in Section |3~21 
we expect only 7 spurious sources to remain in our final source 
catalogue. Measurements of the fraction of flux amplified sources 
in our simulations indicate that confusion is unimportant (< 2% of 
our sources are confused in all energy bands) at the fluxes probed 
by this survey. This is confirm ed via our a cro ss-correlation with 
the Chandra catalogue of McHardv et al. (2003). We find that only 
2% of our sources have ambiguous Chandra counterparts. 

The average input and output fluxes from our simulations at 
bright fluxes agree to within 1%. This has important consequences 
for the determination of source counts, and suggests any systematic 
deviation between input and output source counts should be very 
small. Our simulations demonstrate that Eddington bias affects our 
measured source counts by at most 23 per cent at the faintest fluxes. 
However, we correct for the effects of Eddington bias when fitting 
our source counts, using the distribution of S out /Sj„ p obtained from 
our simulations, as described in Section l3~3l 

Our simulations indicate that the offset in position between 
input and output sources is less than 2" in 68% of cases within 
the central 9' of the XMM-Newton field of view fSection l3~TY We 
therefore expect the majority of our sources to have positions ac- 
curate to 2". However, the EMLDETECT positional errors output 
from the detection chain are too small and an additional systematic 
positional error of 0.76" is required (as determined in Section l43l 
to match the real and simulated positional error distributions. The 
high fidelity of our source positions and fluxes, coupled with in- 
significant confusion illustrate the high reliability of the survey. We 
are able to construct useful X-ray spectra down to a 2-5 keV band 
flux of 3 x 1(T 15 erg citT 2 s , well beyond the initial survey goals 
IPaee et alferjol . 

However, in order to extract the maximum scientific informa- 
tion from the survey, the X-ray sources need to be optically identi- 
fied, and compared with sources detected in other wavebands, such 
as our radio source catalogue ISevmour et all2004l) . This requires 
highly accurate source positions to identify the correct counter- 
part. The additional systematic positional error has important con- 
sequences for optical identification. Using the optical source counts 



determined directly within the 13 H field i McH ardv etall2003h wc 
have calculated the expected number of chance coincidences within 
an appropriate radius of each source within the field (corresponding 
to the 95% positional error radius of each source). Table shows 
the expected number of chance coincidences for the 225 XMM- 
Newton sources as a function of optical magnitude; it also shows 
the actual number of optical sources observed within the error cir- 
cles. X-ray/optical associations at magnitudes R < 23 are firm, but 
approximately half of those at magnitudes 23 < R < 24 are likely to 
be spurious. At fainter magnitudes the majority of optical counter- 
parts are probably unrelated to the X-ray source. This highlights the 
importance of the Chandra coverage over the field which provides 
accurate X-ray positions and extends the range in which we can 
correctly select an optical counterpart to a magnitude of R = 26. 
This is an issue that must be considered in any deep XMM-Newton 
survey since if ignored it may lead to a bias in optical identifica- 
tions. Chandra source counts determined by source type find an 
increasing contribution of absorbed AGN and normal galaxies at 
the fainter fluxes (Bauer et al. 2004). This suggests that broad line 
AGN, which dominate at brighter fluxes, are most likely to be cor- 
rectly identified even without Chandra positions. However, in an 
XMM-Newton survey of similar or greater depth than ours, Chan- 
dra coverage is essential to identify correctly the optical counter- 
parts of the faintest sources, which are more likely to be associated 
with absorbed AGN or normal galaxies. 

We have carr ied out a similar com parison with the 13 w field 
radio catalogue of lSevmour et alJ 120041) . The predicted number of 
chance coincidences (as predicted from the best-fit starburst and 
AGN population model of lSevmour et all2004h are compared with 
the observed number of radio counterparts in Table [8] We are able 
to securely identify the correct radio counterpart to a radio flux of 
O.lmJy without Chandra positions. Below this flux, we need Chan- 
dra positions in order to reliably identify real associations between 
X-ray and radio sources. The relative ease of the identification of 
radio counterparts compared with that of optical counterparts is due 
to the paucity of the radio source population. 

Eddington bias most affects our counts at faint fluxes, where 
the statistical errors on our fluxes are largest. For XMM-Newton 
surveys with shorter exposure times, source counts at faint fluxes 
will be significantly affected by Eddington bias. Eddington bias 
at a particular flux depends on the relative errors on the flux, 
and upon the slope and n ormalisation of the source counts dis- 
tribution l Eddin gtor][T913l: lTeerikorpill2004 . At a flux of 10~ 15 
erg cirr 2 s" 1 the relative difference in signal to noise between ex- 
posures of 120 and 40 ks is ~ 1.7. Given that our source counts are 
boosted by ~ 10% at 3 x 10~ 15 erg cirr 2 s we expect the source 
counts in a survey with a 40 ks exposure time to be boosted by 
~ 29% at 3 x 10~ 15 erg cirr 2 s _I . This is of particular relevance to 
shallo wer surveys such as the HELLAS2XMM survey iBal di et all 
2002), which should be heavily affected by Eddington bias at faint 
fluxes. 

5.3 Contribution to the XRB 

In this section we examine what fraction of the XRB we can probe 
with our survey by comparing the integrated emission from our 
source counts with various measurements of the XRB intensity. 
To obtain the total 1-2 keV emission from point sources, we in- 
tegrate our 0.5-2 keV differential source counts between fluxes of 
2 x 10~ 16 erg cirr 2 s~' and 3 x 10~' 4 erg cirr 2 s _I ; for fluxes greater 
than 3 x 10~ 14 erg cnr 2 s -1 we h ave add ed the integrated counts ob- 
tained from the RIXOS survey <Mason et alj2000l) . In the 1-2 keV 
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Figure 14. Discrete source contribution to the 1-2 keV X-ray background 
from 0.1 - 1.0 X 10~' 3 erg cirT 2 s~'. The solid curve shows the integrated 
emission from the best fit to our source counts, the dotted curves indicate 
the lower and upper 95% error on the integrated emission. Measurements 
of the XRB intensity taken from the literature are overlaid and labelled for 
comparison. Our source counts account for 47 - 92% of the XRB intensity 
in this energy band. 



Figure 15. Discrete source contribution to the 2-5 keV X-ray background 
from 0.1 - 1.0 X 10~' 3 erg cirT 2 s~'. The solid curve shows the integrated 
emission from the best fit to our source counts, the dotted curves indicate 
the lower and upper 95% error on the integrated emission. Measurements 
of the XRB intensity taken from the literature are overlaid and labelled for 
comparison. Our source counts account for 42 - 92% of the XRB intensity 
in this energy band. 



band our source counts account for 51-100% of the XRB intensity. 
Our contribution to the XRB as a function of flux is illustrated in 
Fig. 1141 The main uncertainty in the resolved fraction lies in the 
range of reported values for the absolute normalisation of the XRB 
which are discrepant by u p to 30%. Th e XRB measurements of 
iGendreau etal] J 19951) and lvecchi etal] J] 999) are taken to repre- 
sent the lower and upper normalisations respectively in this energy 
band. 

In Fig.EU we show the contribution of our sources to the inte- 
grated 2-5 keV XRB intensity. We integrated our differential source 
counts between fluxes of 2 x 10~ 15 erg citT 2 s _i and 1 x 10" 13 
er g cm' 2 s" 1 but this time we added on t he integrated counts 
of ICagnoni. Delia Ceca & M accacarol ll998l). With r espect to the 
XRB measurements of lGeTdreau^^lT i ^95]) and IVecchi et alJ 
( 1999), our source counts contribute 50-93% of the XRB inten- 
sity in the 2-5 keV band. Our source counts would account for the 
whole 2-5 keV XRB intensity measured bv lMarshall et alJ <198fjl) . 
Therefore, the XRB intensi ty within the 13 " field must be higher 
than that measured by Mar shall et alj fT980). Indeed, given that our 
source counts are rich around S k„ ee , the XRB intensity within the 
13" field may be somewhat higher than the global average. 

We have high quality X-ra y spectra to 2-5 keV fluxes of 3 x 
10~ 15 erg cirr 2 s~' (Page et al. 2003, 2005). At these fluxes we have 
already resolved ~ 30% of the XRB, and hence the 13" survey is 
readily suitable for the study of the nature and physical properties 
of the major XRB producing populations. 



6 CONCLUSIONS 

In this paper we have presented the complete catalogue of 225 
sources in the XMM-Newton 13" deep field which covers a sky 



area of 0.18 deg 2 centred on position 13h 34m +37° 53'. We reach 
source densities of 700 deg~ 2 at a flux of 4. 1 x 10~ 16 erg cirr 2 s in 
our lowest energy band. In the 0.5-2 keV band we find source den- 
sities of 1300 at a limiting flux of 4.5 x 10~ 16 erg enr 2 s -1 . At 
harder energies we reach fluxes a factor ~ lOx brighter, with source 
densities of 900 and 300 deg~ 2 at limiting fluxes of 1.1 x 10~ 15 
erg enr 2 s and 5.3 x 10~ 15 erg cirr 2 s in the 2-5 keV and 5-10 
keV energy bands respectively. 

We have carried out extensive simulations of the detection pro- 
cess in order to assess the reliability of our source catalogue. Sim- 
ulations indicate that confusion is small (< 2%) in our images. We 
have curtailed the sourcelist directly derived from the SAS accord- 
ing to the optimal statistical detection likelihoods in each band de- 
termined from our simulations. We expect < 7 spurious sources to 
remain in the final catalogue. 

Within the central 9' of the XMM-Newton field of view, posi- 
tional errors are less than 2" for 68% of our simulated sources and 
our input/output fluxes agree to within 1% at bright fluxes. Com- 
parison of the input/output source positional offsets from our simu- 
lations with the positional offsets found between the XMM-Newton 
and Chandra counterparts suggest that an additional systematic 
error of 0.76" should be added in quadrature to the EMLDE- 
TECT positional errors. This poses no problem for radio counter- 
pa rt id entification. We match our catalogue with the radio catalogue 
of Sevmour et al. 1 2004) and are confident in our choice of radio 
counterpart to a flux of O.lmJy. In the optical we are confident in 
our choice of optical counterparts to a magnitude of R = 23 us- 
ing our XMM-Newton positions. For magnitudes fainter than this 
we need Chandra positions w hich allow reliable id entification of 
optical counterparts to R = 26 iMcHardv et al. 2003). 

We have computed the best-fit parameters for the differential 
N(S ) function using a method adapted from [Murdoch et aT]ll973l) . 
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In all but our hardest energy band the data are described by a dou- 
ble powerlaw model with a break flux at ~ 10~ 14 erg cnT 2 s _1 . 
The counts below the break are flatter than the Euclidean case. 
Our measured source counts are in broad agreemen t with pre- 
vious determinations from the Chandra deep fields I Brandt et al. 
| 200 it IRosati et alj I2002I) . XMM-Ne wton iHasinger et alJ 12001b 
baldi et alJ 120021) and BeppoSAX iFiore et alJ l200ll) surveys. 
The overall normalisation in the 0.5-2 keV band is similar to, 
thoug h higher than, that fo und in the CDF-N survey iBrandt et alJ 
l200ll). Both the CD F-S iRosati et alJl2002l) and Lockman hole 
( Hasinaer et al. 2001) normalisations are significantly lower than 
found here. This field-to-field variation in source density may be 
attributed to the underlying clustered large scale structure which 
imprints a variation in source density on small scales. 

In the 2-5 keV ba nd our faint end normalisation is consis- 
tent with the CDF-N iBrandt et all 1200 ll) . CDF- S iRosati et al 
I2OO2I) and XMM-Newton Lockman hole counts iHasinger et al 
200 1 ) . There is minimal flux overlap with the ASCA su rveys 
ICagnoni. Delia Ceca & Maccacaro! 1998blUeda et alJl999lal) . The 
2-5 keV band ASCA counts lie above tho se found in this survey 
ICagnoni. Delia Ceca & Maccacaro! fl9 98) although they are still 
consistent with our model fits. In our hardest energy band (5-10 
keV) our counts are in broad agreement with those from BeppoSAX 
and ASCA , and with the 5-10 keV counts in the Lockman hole 
( Hasin ger et alj200ll) . Again there is little overlap with the brighter 
ASCA and BeppoSAX surveys. 

The sources in our survey straddle the break in the source 
counts. Accounting for the uncertainty in our source counts and the 
absolute normalisation of the XRB we resolve 51-100(50-93)% of 
the 1-2 (2-5) keV XRB emission. At the break in the source counts 
we resolve ~ 30% of the 2-5 keV XRB. At these fluxes we have 
X-ray spectra and are therefore able to study the emission mech- 
anisms and investigate absorption in a significant fraction of the 
X-ray source population. 
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Table 9: 13 H XMM-Newton deep field catalogue. All fluxes are in units of 10 
erg cm -2 s _I . The columns labelled 'Counts' and 'Counts error' give the best fit num- 
ber of source counts and the lcr error as output from EMLDETECT. Flux errors are at lc 
confidence. If in any energy band there are fewer counts at the source position in the im- 
age than in the background map, then the number of source counts and the source flux are 
recorded as zero for that band. In this case the 'Counts error' and 'Flux error' columns rep- 
resent the lcr upper limit. The 'Pos error' column gives the EMLDETECT positional error 
(lcr) in arc seconds. The 'Offaxis angle' column gives the XMM-Newton offaxis angle in 
arc minutes. The 'Cha ndra No.' column gives the Chandra catalogue number taken from 
iMcHardv etal] l2003l) . The 'Chandra offset' column gives the offset in position between 
the XMM-Newton and Chandra source, in arcsec. 
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13.21 


5.53 


46 


2.77 


18 


13 33 43.79 


37 54 55.76 


2.56 


9.94 


12.00 


6.49 


34.76 


9.48 


0.71 


4.00 


2.67 


5.62 


0.30 


0.16 


0.81 


0.22 


0.04 


0.21 


0.59 


1.25 


1.74 


1.29 






19 


13 33 44.43 


37 57 53.49 


0.38 


10.26 


353.57 


21.39 


422.51 


23.96 


103.76 


13.33 


11.93 


8.58 


9.03 


0.55 


9.90 


0.56 


5.62 


0.72 


2.69 


1.94 


27.24 


2.21 


29 


1.95 


20 


13 33 45.45 


37 58 08.04 


0.92 


10.14 


128.52 


15.44 


114.80 


15.21 


14.87 


8.98 


9.44 


8.43 


3.25 


0.39 


2.66 


0.35 


0.80 


0.48 


2.11 


1.88 


8.82 


2.01 






21 


13 33 46.61 


38 00 21.56 


0.74 


10.85 


21.83 


7.18 


170.94 


15.57 


62.32 


11.11 


11.40 


7.61 


0.69 


0.23 


4.89 


0.45 


4.06 


0.72 


3.16 


2.11 


12.80 


2.28 


78 


0.76 


22 


13 33 47.52 


37 53 51.50 


1.11 


9.27 


33.80 


8.58 


45.73 


10.90 


0.00 


7.31 


8.14 


8.38 


0.80 


0.20 


1.00 


0.24 


0.00 


0.37 


1.66 


1.71 


3.46 


1.78 







23 


13 33 48.28 


37 53 33.71 


0.64 


9.16 


13.48 


6.96 


50.72 


10.11 


199.75 


16.71 


48.63 


10.97 


0.32 


0.17 


1.11 


0.22 


10.05 


0.84 


9.97 


2.25 


21.45 


2.42 


76 


0.66 


24 


13 33 48.60 


37 58 07.79 


0.57 


9.55 


146.08 


14.13 


166.80 


15.63 


43.17 


9.69 


5.77 


6.80 


3.56 


0.34 


3.75 


0.35 


2.25 


0.50 


1.22 


1.44 


10.78 


1.61 


40 


1.25 


25 


13 33 50.68 


37 49 46.61 


0.85 


10.01 


19.02 


7.54 


132.61 


14.58 


45.80 


10.49 


3.45 


5.66 


0.48 


0.19 


3.05 


0.34 


2.44 


0.56 


0.75 


1.24 


6.72 


1.41 


128 


1.36 


26 


13 33 50.80 


37 57 16.61 


1.77 


8.88 


8.66 


6.42 


25.68 


8.03 


32.51 


9.23 


11.06 


7.91 


0.20 


0.15 


0.53 


0.17 


1.56 


0.44 


2.15 


1.53 


4.43 


1.61 






27 


13 33 53.59 


38 02 03.64 


0.77 


10.73 


42.68 


8.30 


117.78 


13.64 


18.24 


7.38 


11.30 


8.55 


1.15 


0.22 


2.91 


0.34 


1.04 


0.42 


2.71 


2.05 


7.81 


2.13 


72 


1.29 


28 


13 33 54.21 


38 02 53.17 


1.23 


11.21 


35.02 


8.06 


16.36 


7.20 


12.39 


6.80 


32.15 


9.36 


1.64 


0.38 


0.67 


0.30 


1.11 


0.61 


13.07 


3.81 


16.49 


3.89 






29 


13 33 54.91 


37 51 25.87 


1.81 


8.50 


17.87 


6.74 


33.18 


9.28 


0.00 


5.35 


0.00 


5.89 


0.39 


0.15 


0.67 


0.19 


0.00 


0.25 


0.00 


1.09 


1.06 


1.14 


175 


3.41 


30 


13 33 55.88 


37 52 57.85 


0.28 


7.81 


227.18 


16.97 


636.62 


26.92 


196.33 


16.74 


75.24 


12.65 


4.87 


0.36 


12.55 


0.53 


8.87 


0.76 


13.42 


2.26 


39.71 


2.46 


17 


1.04 


31 


13 33 57.90 


37 49 59.01 


2.91 


8.70 


20.13 


7.65 


25.94 


9.21 


0.00 


2.37 


0.00 


2.55 


0.45 


0.17 


0.53 


0.19 


0.00 


0.11 


0.00 


0.48 


0.98 


0.56 






32 


13 33 58.56 


37 59 37.55 


0.23 


8.46 


478.40 


23.67 


1031.89 


34.31 


320.33 


20.81 


70.33 


12.92 


10.84 


0.54 


21.58 


0.72 


15.50 


1.01 


13.59 


2.50 


61.51 


2.84 


7 


0.75 


33 


13 33 58.87 


38 00 24.88 


1.17 


8.87 


13.57 


6.81 


30.94 


8.16 


22.51 


7.78 


14.47 


8.31 


0.44 


0.22 


0.90 


0.24 


1.43 


0.50 


3.83 


2.20 


6.60 


2.28 






34 


13 33 59.21 


38 05 58.32 


1.37 


13.03 


7.77 


4.94 


17.99 


5.00 


18.54 


6.00 


1.12 


3.23 


1.06 


0.67 


1.88 


0.52 


3.62 


1.17 


1.21 


3.49 


7.78 


3.78 
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3.20 
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0.2-0.5 keV 


0.5-2.0 keV 


2.0-5.0 keV 


5.0-10.0 keV 


0.2-0.5 keV 


0.5-2.0 keV 


2.0-5.0 keV 


5.0-10.0 keV 


0.2-10.0 keV 
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35 


13 33 59 56 

± J JJ y . ^>\j 


37 51 39.75 


2.34 


7.57 


5.28 


6.21 


36.92 


10.18 


9.90 


7.53 


5.15 


6.68 


0.11 


0.12 


0.68 


0.19 


0.42 


0.32 


0.86 


1.11 


2.06 


1.18 






36 


13 33 59.63 


37 49 26.35 


12.19 


8.75 


0.00 


4.59 


90.42 


19.99 


0.00 


5.21 


0.00 


7.12 


0.00 


0.10 


1.86 


0.41 


0.00 


0.25 


0.00 


1.35 


1.86 


1.44 






37 


13 34 00.04 


37 49 12.27 


0.32 


8.84 


373.03 


20.83 


712.58 


28.63 


131.11 


14.00 


3.02 


5.26 


8.45 


0.47 


14.86 


0.60 


6.31 


0.67 


0.58 


1.01 


30.21 


1.44 


4 


1.49 


38 


1 3 34 00 9 1 

±j J7 \j\j.y± 


38 01 25 24 


0.95 


9.24 


5.46 


5.57 


88.41 


12.16 


44.95 


10.00 


8.60 


7.99 


0.13 


0.13 


1.92 


0.26 


2.25 


0.50 


1.75 


1.63 


6.06 


1.73 


84 


0.51 


39 


13 34 01 05 

±J J7 \J ± . WJ 


37 54 03 99 

j / J 7 \j*j.yy 


0.35 


6.60 


244.79 


17.97 


549.22 


25.28 


261.35 


19.83 


76.08 


13.84 


4.95 


0.36 


10.35 


0.48 


11.50 


0.87 


12.64 


2.30 


39.45 


2.53 


23 


1.00 


40 


13 34 01 06 


38 00 25 75 

jo ynj z. j . / j 


1.50 


8.55 


6.96 


6.07 


19.52 


7.67 


58.74 


10.85 


11.63 


7.96 


0.17 


0.15 


0.43 


0.17 


2.94 


0.54 


2.38 


1.63 


5.92 


1.73 






41 


13 34 01 26 


37 46 47 56 

J 1 " / . JU 


0.90 


10.41 


8.67 


6.83 


147.83 


15.54 


59.26 


11.84 


14.80 


8.84 


0.23 


0.18 


3.55 


0.37 


3.31 


0.66 


3.40 


2.03 


10.48 


2.18 


48 


1.05 


42 


1 3 34 02 62 

± J J*T VJZ..UZ. 


37 51 29 14 


0.36 


7.12 


144.04 


14.09 


419.12 


22.19 


99.91 


13.25 


21.77 


9.24 


2.80 


0.27 


7.52 


0.40 


4.12 


0.55 


3.47 


1.48 


17.92 


1.65 


43 


0.88 


43 


13 34 03 09 

±J J7 \J^J.\Js 


37 53 21 42 


0.50 


6.34 


2.26 


4.13 


95.74 


12.95 


287.86 


19.70 


90.74 


13.78 


0.04 


0.08 


1.62 


0.22 


11.20 


0.77 


13.50 


2.05 


26.36 


2.20 


52 


0.58 


44 


1 3 34 07 1 S 

± J J7 U/.IO 


37 59 59 34 


1.23 


7.36 


18.32 


10.34 


0.00 


3.70 


45.94 


14.61 


15.40 


10.88 


0.37 


0.21 


0.00 


0.07 


1.93 


0.62 


2.57 


1.81 


4.87 


1.93 






45 


13 34 08.05 


38 06 26.62 


0.92 


12.62 


33.53 


7.81 


94.69 


12.52 


6.79 


7.46 


33.96 


10.12 


1.29 


0.30 


3.29 


0.43 


0.54 


0.59 


12.00 


3.58 


17.12 


3.67 


60 


1.39 


46 


1 3 34 08 36 


37 52 19 21 

J / J _ 1 .:/ . ^ 1 


0.50 


5.72 


115.74 


12.85 


256.77 


18.53 


73.97 


11.86 


28.75 


9.60 


2.01 


0.22 


4.12 


0.30 


2.71 


0.43 


4.00 


1.34 


12.84 


1.45 


68 


1.08 


47 


1 3 34 08 39 

± J J7 V /O . J y 


37 47 50 85 


2.69 


8.72 


5.54 


6.28 


15.49 


8.33 


34.65 


10.02 


13.50 


8.48 


0.12 


0.14 


0.32 


0.17 


1.66 


0.48 


2.56 


1.61 


4.67 


1.69 






48 


13 34 08.59 


37 54 40.76 


0.53 


5.06 


94.30 


12.27 


234.91 


18.41 


118.64 


15.53 


13.93 


8.44 


1.56 


0.20 


3.58 


0.28 


4.10 


0.54 


1.83 


1.11 


11.08 


1.28 


58 


1.38 


49 


1 3 34 08 72 

± J J7 V/O. / 


38 03 49 23 

jo v/j T\y . ^_ j 


0.47 


10.22 


112.09 


12.97 


331.70 


21.03 


43.75 


10.19 


12.05 


9.31 


2.94 


0.34 


7.99 


0.51 


2.44 


0.57 


2.80 


2.16 


16.16 


2.32 


25 


1.01 


50 


1 3 34 08 79 

1 J J*T I/O. / 7 


37 57 06 43 


0.27 


5.47 


249.06 


17.77 


860.91 


31.47 


287.71 


20.38 


88.02 


13.60 


4.29 


0.31 


13.63 


0.50 


10.37 


0.73 


12.13 


1.87 


40.42 


2.10 


28 


0.45 


51 


1 3 34 09 96 


37 54 3 1 62 


0.78 


4.80 


0.82 


3.39 


23.43 


9.06 


128.60 


14.47 


47.77 


11.24 


0.01 


0.06 


0.35 


0.14 


4.38 


0.49 


6.15 


1.45 


10.89 


1.53 


141 


0.82 


52 


1 3 34 10 61 


37 59 56 03 

J 1 J7 JU.UJ 


0.31 


6.84 


641.48 


29.15 


351.57 


22.04 


94.69 


12.23 


15.53 


9.25 


17.53 


0.80 


8.66 


0.54 


5.18 


0.67 


3.43 


2.04 


34.80 


2.35 


9 


0.46 


53 


13 34 11.42 


37 47 57.32 


0.51 


8.29 


172.60 


15.68 


234.29 


18.23 


62.26 


10.90 


20.36 


9.54 


3.72 


0.34 


4.66 


0.36 


2.85 


0.50 


3.68 


1.73 


14.92 


1.87 


66 


1.76 


54 


1 3 34 1 2 36 

±J J7 ± i. . J v/ 


37 59 1 o 36 


2.01 


6.06 


28.74 


11.66 


0.00 


3.54 


0.00 


5.63 


3.76 


5.90 


0.53 


0.22 


0.00 


0.06 


0.00 


0.22 


0.56 


0.88 


1.10 


0.94 






55 


13 34 13.03 


37 58 31.63 


0.77 


5.52 


81.50 


11.53 


147.43 


14.16 


43.77 


9.56 


10.94 


8.25 


1.85 


0.26 


3.02 


0.29 


1.97 


0.43 


1.94 


1.47 


8.79 


1.58 


85 


0.71 


56 


13 34 13.59 


37 45 39.06 


1.22 


10.12 


32.55 


8.52 


65.93 


10.90 


0.00 


2.05 


0.00 


5.19 


0.84 


0.22 


1.57 


0.26 


0.00 


0.11 


0.00 


1.18 


2.42 


1.23 






57 


1 3 34 1 3 62 

±J J"T ± J . U 


38 07 23.60 


0.89 


13.12 


23.57 


7.54 


195.39 


17.42 


69.29 


11.80 


13.87 


9.29 


0.83 


0.26 


6.28 


0.56 


5.18 


0.88 


4.59 


3.08 


16.89 


3.26 


88 


1.66 


58 


13 34 14.15 


38 04 41 20 


1.00 


10.54 


9.96 


6.36 


75.06 


12.90 


95.57 


13.20 


29.41 


9.64 


0.27 


0.17 


1.86 


0.32 


5.46 


0.75 


7.09 


2.32 


14.67 


2.47 


132 


1.35 


59 


13 34 14 37 

± J J7 ± 7 . J / 


37 52 32 09 

J / J^. -J ^- .\J y 


1.43 


4.59 


0.00 


1.14 


37.94 


8.56 


19.46 


8.35 


0.00 


6.75 


0.00 


0.03 


0.84 


0.19 


0.93 


0.40 


0.00 


1.27 


1.77 


1.34 


119 


1.03 


60 


13 34 14 79 

1 J 17. / 7 


37 59 59 67 


0.85 


6.36 


50.74 


10.19 


130.43 


14.39 


22.30 


8.86 


17.50 


8.82 


0.97 


0.19 


2.29 


0.25 


0.89 


0.35 


2.71 


1.37 


6.86 


1.45 


120 


1.20 


61 


13 34 14.81 


37 51 30 81 

J / J 1 JVJ.ul 


0.38 


5.12 


206.64 


15.92 


353.49 


21.06 


0.00 


2.77 


0.00 


3.77 


3.48 


0.27 


5.49 


0.33 


0.00 


0.10 


0.00 


0.50 


8.98 


0.67 


81 


0.12 


62 


1 3 34 I 5 49 


38 03 05 78 

JO V/ J V/ J . / O 


1.80 


8.97 


11.63 


6.18 


13.12 


7.94 


32.35 


9.58 


13.20 


8.16 


0.27 


0.14 


0.28 


0.17 


1.60 


0.47 


2.66 


1.64 


4.81 


1.72 






63 


13 34 15.54 


37 52 27.80 


1.31 


4.43 


0.00 


1.52 


43.88 


14.26 


65.88 


11.49 


19.04 


9.60 


0.00 


0.03 


0.68 


0.22 


2.32 


0.40 


2.53 


1.28 


5.54 


1.36 






64 


13 34 17.01 


37 59 48 75 

J / J^/ 'TO. / J 


0.96 


5.95 


48.10 


9.74 


87.75 


12.41 


22.14 


8.52 


1.93 


5.23 


0.87 


0.18 


1.46 


0.21 


0.84 


0.32 


0.28 


0.77 


3.45 


0.87 


170 


0.72 


65 


13 34 17.51 


37 57 21.91 


0.14 


4.10 


1590.56 


41.23 


2266.48 


49.01 


681.39 


28.69 


136.27 


15.72 


25.25 


0.65 


33.02 


0.71 


22.45 


0.95 


16.96 


1.96 


97.67 


2.38 


5 


0.64 


66 


13 34 1 7 SO 


37 44 31.39 


2.04 


10.89 


3.97 


5.25 


23.78 


9.25 


43.92 


10.78 


0.00 


7.79 


0.12 


0.15 


0.64 


0.25 


2.76 


0.68 


0.00 


2.03 


3.51 


2.16 






67 


1 3 34 1 8 87 

± J ±0.0/ 


37 58 56 69 

j i jo *->\j .\j y 


0.32 


5.04 


18.29 


6.88 


403.89 


22.27 


415.62 


22.78 


110.89 


13.89 


0.31 


0.12 


6.34 


0.35 


14.87 


0.81 


15.03 


1.88 


36.56 


2.08 


63 


0.76 


68 


13 34 19.24 


37 43 49 99 

J 1 7 J ^y.yy 


1.55 


11.47 


0.00 


3.55 


37.29 


9.86 


32.42 


9.24 


17.94 


9.22 


0.00 


0.11 


1.02 


0.27 


2.08 


0.59 


4.84 


2.49 


7.94 


2.57 


137 


1.79 


69 


1 3 34 1 9 26 

±J J ■ ± V . i. V/ 


37 50 29 60 


0.52 


5.32 


105.10 


13.16 


217.42 


17.46 


67.50 


11.43 


19.94 


9.27 


1.87 


0.23 


3.56 


0.29 


2.51 


0.43 


2.81 


1.31 


10.75 


1.42 


116 


1.40 


70 


13 34 19 56 

±J ±_/.J\_l 


37 51 47 77 

J/ J ± *T / « / / 


1.40 


4.26 


0.00 


5.64 


22.75 


8.79 


43.62 


10.11 


14.49 


8.37 


0.00 


0.09 


0.33 


0.13 


1.45 


0.34 


1.80 


1.04 


3.58 


1.10 






71 


13 34 19.97 


37 54 00.29 


0.61 


2.96 


23.51 


7.43 


170.55 


15.77 


63.23 


10.82 


27.08 


8.46 


0.41 


0.13 


2.70 


0.25 


2.23 


0.38 


3.60 


1.12 


8.94 


1.22 


93 


0.66 


72 


13 34 20.80 


37 59 30.74 


0.96 


5.30 


11.54 


6.63 


25.22 


8.81 


91.13 


12.72 


24.21 


8.82 


0.20 


0.12 


0.41 


0.14 


3.38 


0.47 


3.41 


1.24 


7.40 


1.34 


191 


1.72 


73 


13 34 20.82 


37 54 59.88 


0.50 


2.64 


99.21 


12.17 


201.23 


16.48 


80.58 


11.46 


11.53 


7.36 


1.59 


0.20 


2.96 


0.24 


2.66 


0.38 


1.43 


0.91 


8.64 


1.03 


59 


1.30 


74 


13 34 9 1 QO 
1 J z 1 . yu 


38 04 50 89 

JO Ut JU.07 


1.34 


10.22 


13.39 


6.83 


70.31 


10.93 


17.24 


8.12 


0.00 


4.95 


0.36 


0.18 


1.71 


0.27 


0.97 


0.46 


0.00 


1.17 


3.03 


1.29 


1 33 


n 61 


75 
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1.73 
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08 
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50 
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U.OO 


1 88 
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196 


0.99 


76 


13 34 22.11 


37 48 03.61 


1.88 


7.26 


6.37 


6.82 


10.70 


7.81 


51.41 


10.64 


29.54 


10.23 


0.13 


0.14 


0.20 


0.14 


2.17 


0.45 


4.84 


1.68 


7.33 


1.75 


214 


1.05 


77 


13 34 22.15 


38 06 20.07 


0.67 


11.66 


74.21 


10.55 


187.68 


15.93 


51.17 


10.51 


5.49 


6.60 
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APPENDIX A: SYNTHETIC FIELD GENERATION 
Al Method 

Input sourcelists were generated independently for each of the four 
energy bands. The differential source count function, N(S) , was 
modelled as a broken power law distribution: 



N(S) 



K\ (S < Sknee) 

^ 2 (sfbr) ^ Sknee) 



(Al) 



The normalisation constants were chosen such that the func- 
tion is continuous at the knee, i.e. K\jKi = (S knee/S normY' 72 ■ 
The parameters selected for the simulations were based on results 
from the CDF-S (see lRosati et all2002t) . A photon index, T = 1,7 
(approximately that of the XRB in the 0.2-10 keV energy range), 
was used to convert from the Chandra energy bands to the energy 
bands used in this study. We have slightly modified the Chandra 
N(S) parameters via an iterative process so that the 'output' source 
counts approximately match those seen in the 13 H field data. 

In order to incorporate the effects of source confusion, we set 
our simulated source flux limits, 5/„„, to values approximately five 
times fainter than the limit reached by the 13 H data. The mean 
number of input sources per field per energy band, (A'), was calcu- 
lated by integrating N(S ) from 5 /,„> to infinity. The actual number 
of sources used in each simulated field was taken randomly from a 
Poisson distribution about (N). For each of these sources, a flux was 
randomly assigned from the appropriate N(S ) distribution. Any in- 
put source having a flux greater than twice that of the brightest 
source in the 13 H data was discarded. This prevented any single 
simulated field being dominated by an extremely bright source: a 
situation not seen in the 13 H data. The effects of source clustering 
are ignored in this analysis, hence each input source was assigned 
a purely random position within the field. Finally, the source fluxes 
were converted to count rates. 
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Figure Bl. Scatter plot of input vs output flux in the 0.5-2 keV energy band. 
We only show those detected sources having DET_ML 5= 5 and valid input 
counterparts within 5", 8", and 10" for input offaxis angles of 0—9', 9-12' 
and > 12' respectively. 



A3 Simulated image analysis 

We used a single-band version of the source searching process (as 
described in Section l2~2l on the simulated images. This process 
was repeated independently in each of the four energy bands. Two 
iterations of the background fitting/source finding process were car- 
ried out per band per field. The output sourcelist was curtailed at a 
detection likelihood of DETJVIL = 5. 



A2 Imaging characteristics 

To convert the simulated input sourcelists to images, one has to 
take account of the complex point spread function (PSF) of the 
EPIC cameras. We used the 'MEDIUM' accuracy PSF description 
( Kirch 2004) from the XMM-Newton calibration files, which is also 
the model used by the SAS source searching task EMLDETECT. 
This PSF description consists of a number of small ray-traced maps 
each describing the PSF at a particular energy and off-axis angle, 
covering the full EPIC field of view and energy range. These maps 
were interpolated in energy and off-axis angle, allowing us to eval- 
uate the fraction of any source's flux within any image pixel. The 
SAS-generated exposure maps from the 13 H field give the effec- 
tive exposure times and vignetting corrections for each energy band 
and EPIC camera. The simulated images were generated pixel by 
pixel by summing the contribution from all input sources. We added 
a two-component, (vignetted and un-vignetted), synthetic back- 
ground to the simulated images to reproduce that observed in the 
13 H field. The correct level of this background was determined 
through an iterative process because undetected, faint simulated 
sources contribute significantly to the diffuse background. Finally, 
the value of each pixel was randomly drawn from the appropriate 
Poisson distribution to simulate photon counting noise. 



APPENDIX B: ASSESSING THE IMPACT OF 
CONFUSION 

Bl Simulation method 

In order to investigate the ultimate confusion limit of deep 
XMM-Newton surveys, we have produced a small number of ex- 
tremely deep simulations, with exposure times 1000X that of the 
13" field. The input sourcelists reach fluxes of a few 10~ 19 - 10" 18 
erg cirr 2 s -1 depending on energy band. No background or Poisson 
noise were added to the images in these simulations, so differences 
between the input and output source counts arise solely as a re- 
sult of confusion. The sources were matched within a cutoff radius, 
r cut , which depends on XMM-Newton offaxis angle as described in 
Section lXTI 

B2 Results 

Fig. IB II shows the distribution of S ml1 /S ,■„,, for our 10 simulations 
in the 0.5-2 keV energy band. The scatter on S ou ,/S inp increases 
dramatically below a flux of 10~ 16 erg cirr 2 s -1 . At similar fluxes, 
the input-output position offsets, shown in Fig. lB2l become much 
larger. The source density at this flux level is approximately 2000 
deg~ 2 and represents the limit beyond which source properties can- 
not be recovered reliably irrespective of exposure time. 
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Figure B2. Greyscale image showing the distribution of positional offsets 
between output and input source locations as a function of input flux in the 
0.5-2 keV energy band. Those sources within the central 9' of the XMM- 
Newton field of view having DET.ML ^ 5 are shown. The concentration of 
positional offsets as a function of Sj np is indicated by the darkness of the 
greyscale image. The three contours show the distances within which 68, 
90, and 95 per of the data lie. 



Fig. IB3I show the average input and output integral source 
counts in each of our four energy bands. The input and out- 
put source counts are already discrepant when the source density 
reaches 2000 deg~ 2 , and diverge rapidly at fainter fluxes. The ap- 
proximate input fluxes corresponding to source densities of 2000 
deg~ 2 are 3 x 10~ 17 , lO" 16 , 2 x 10" 16 and 3 x lO" 16 erg cnr 2 s _1 in 
the 0.2-0.5, 0.5-2, 2-5 and 5-10 keV energy bands respectively. 
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Figure B3. Comparison of input (dotted line) and output (solid line) simulated source counts. Input source counts were simulated to fluxes of 5 X 10" 
erg cirT 2 s' 1 in the 0.2-0.5 keV and 0.5-2 keV energy bands. In the 2-5 keV and 5-10 keV energy bands the input counts were simulated to fluxes of 1 X 10" 
and 5 X 10~ 18 erg era" 2 s" 1 . Confusion is evident at faint fluxes where the input and output source counts rapidly diverge. 
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